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ABSTRACT 


This thesis is concerned with an investigation 
of the par one OfPehe various structure, factors, 
especially of their long wavelength (wave number q> 0) 
limits, for different types of mixtures. 

First we examine the behaviour of the q-> 0 


limits of the concentration fluctuations S (qi 0 lecand 


Ce 
the various partial structure factors aya Oh ising the 
model of conformal solutions. The model is applicable 
to systems in which (a) the two types of atoms (A and 
B) are of roughly the same size, and (b) the difference 
Ww in pairwise interaction energies between the atoms is 
small (for an ideal agin doe We 0). ties founds that 
the interchange energy w is directly related to Sac (0) 
or the combination a, , (0) + Ay_ (0) - 2a,5 (0). The model 
is applied to calculate an5 0) for Na-K alloys and good 
agreement with experiment is found. The behaviour of 
Sac (0) for alloys where the assumption (a) above does 
not hold is also discussed, using Flory's approximation. 
Next, a phenomenological model is proposed to 
calculate Sac (0) and other thermodynamic properties of 
systems in which the interatomic interactions are strong 
and which form compounds in the solid phase at one or 
more well defined stoichiometric compositions specified 


by A B,) (i,v small) integers). It isgassumed that the 
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binary liquid mixture consists of A and B atoms and 
their chemical associations AB, in chemical equili- 
brium. The formulation is given in general terms, but 
the actual calculations are made by assuming that only 
one type of chemical associations is formed. By 
assuming first, for simplicity, that) the mixture, of 
A band A By can be considered ideal, it is shown that 
the concentration dependence of Sac (0) and hence of the 
various oe depend characteristically on the values 
Gta(tgv). “The iresults)tor Sac (0) Lor, (iy a= (0S 7203, 
(3,1) and (2,1) are in qualitative agreement with those 
determined from measured thermodynamic activity data 
for Mg-Bi, Ag-Al and Hg-K. Itis then shown by numerical 
calculations for four systems(Mg-Bi, Tl-Te, Ag-Al and 
Cu-Sn) that quantitative agreement with experiment can 
be achieved if the ternary mixture of A, B and A, B,, is 
treated in either the conformal or Flory's approximations, 
the interaction parameters being determined from the 
observed free energy. It is concluded that evidence 
for the formation of chemical complexes should be sought 
by making neutron scattering experiments. 

Finally, a simple modification (based on the 
random phase approximation) of the usual hard sphere 
calculations of the structure factors, at arbitrary q, 


ise proposed) to takewinto account the ‘effect: of Lhe wong 


range attractive interactions in weakly interacting 
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systems. The method is first applied to calculate 

S(q) for liquid argon obtaining good agreement with 
experiment. Next, some numerical examples are given 

for mixtures. One finds that the effect of the 
attractive interactions may be expected to be limited 

to q values in the range 0<q <q: where qo is the 
position of the main peak in aj4 (a). The significance 
of these results to the calculation of electrical resis- 


Pivity,1ecc a usuorierfly discussed. 
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The concentration fluctuations Scc as 
a function of concentration ina 
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Ciations of the AB type (y=v=1). 


eae as a function of concentration c 
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Scc as a function of concentration c 
of A atoms in a binary mixture forming 
chemical associations of the AB, type 
(u=3, v=1). 


Sc@ras agkunctioniof concentration c 
of A atoms in a binary mixture forming 
chemical associations of the AB. type 
(N=25.4N=))40 ans 


he parcial ystructunevilactors) ao. asa 
Eunctionvof.concentration (ci-of Ay atoms 
in a binary mixture forming chemical 
associations of the AB type. 


The partial structure factors ajj asa 
function. of :concentration cof Aratoms 
in a binary mixture forming chemical 
associations of the type A.B, CH=37 
v=2). 


The partials structure factors) aj. as a 
function of concentration c of A atoms 
in a binary mixture forming chemical 
associations of the type A,B, (u=2, 
v=1). 


Gy/V RT and n3/V as functions of con- 
Centratizonetor Ti-Te. 


Gu/V RT and n3/N as functions of. con- 
centration for Mg-Bi. 


Gy/VYRT and n3/V” as functions of con- 
centration for Ag-Al. 


Gy/VRT and n3/W as functions O21 (con— 
centrations LoreCu-on. 
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6.3a 


6.4b 


The logarithm of the activity, ln ar. 
versus concentration for Tl-Te. 


Me tlogarrenm of the activity, in 
versus concentration for Mg-Bi. 


the Vogarithm of the activity, in Aangr 
versus concentration for Ag-Al. = 


whe rogarithm of the-activity, In acy: 
versus concentration for Cu-Sn. 


The concentration fluctuation Scc 
versus concentration for Tl-Te. 


The concentration fluctuation Scc 
versus concentration for Mg-Bi. 


The concentration fluctuation Scc 
versus concentration for Ag-Al. 


The concentration fluctuation Scc 
versus concentration for Cu-Sn. 


Heat and entropy of mixing for T1l-Te 
system. 


Sketch of the model potential showing 
d(r) of Lennard-Jones form and x(r) 
as well as the hard core. 


Model S(q) for ¢€/kp = 123°K, 

Y= 93 4x1 OS cempand 40-00-86 er. 
compared with experimental S(q) for 
Arsat 85°h. 


Model 6(q), Scc(q) and A(q) for two 
concentrations (c = 1/2, c= 17/5) . 
taking (Ri=)0010,0T = 65° GC .es ergtical 
parameters). 


Model a;~(q) for the same parameters 
as Fig. (62ca- 

Model 9(q), Scc(q) and A(q) at an © 
elevated temperature: T = 3/2 Tg with 
Ree Oe Sack). 


Model a;.(q) for the same parameters 
as Fig. Geta. 
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Figure 


6.5a 


Model 8(q), Scc(q) and A(q) for 
R = (=-0.2, T = 85°K for two con=- 
centrations. 


Model a,;.(q) for the same parameters 
aseFig, “6 75a. 


Seeiq) tor, two values of R.at tT =—¢85 "5. 
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CHAPTER 1 


INTRODUCTION 


This thesis is concerned with a study.of, 
firstly, the long wavelength (wavenumber q~>0) limits 
of the structure factors of binary molten alloys whose 
constituent atoms interact with one another either 
weakly or strong, and, secondly, the structure factors 
at all q for weakly interacting systems. 

The structure factors play an important role 
in our understanding of interatomic interactions, on 
the one hand, and of various transport properties of 
molten alloys, on the other; for an eee ave review, 
see Faberr41972) 

The structure of (relative arrangement of 
atoms in) a pure liquid is described by a single 
structure factor S(q) which is essentially the Fourier 
transform of the pair distribution function g(r) 

[see Chapter 2]. Since X-ray or neutron scattered 
Sntenscities aze»directly proportional. to S(q), its 
experimental determination is a straightforward matter. 
On the theoretical side, there have been many attempts 
to calculate S(q); notable are those using the hard 
sphere model in the Percus-Yevick approximation 
[Ashcroft and Lekner (1966), etc.] and modifications 


thereof [Woodhead-Galloway and Gaskell (1968), etc.]. 
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For a comprehensive list of references see, for example, 
Croxton (1974). In the zero-wavenumber limit, however, 
S(0) can be exactly expressed in terms of certain 
thermodynamic quantities via the Einstein-Smoluchowski 


or Ornstein-Zernike relation: 
N 


where Kn is the isothermal compressibility, N is the 
total number of atoms in volume V, T is the temperature, 
and kp is the Boltzmann constant. 

In a binary A-B mixture there are three inde- 
pendent structure factors corresponding, basically, to 
the correlations between A-A atom pairs, B-B atom pairs 
and A-B atom pairs. There are several equivalent 
(linearly related) sets of structure factors defined 
in the literature: for example, the partial structure 
factors aj4‘d), or 554 a) [Faber and Ziman (1965), 
Asncrott,. and, Langretn (1967)) ,~and. the number—concen- 
tration-structure factors Syn (4) + Sue (@) and Sac (q) 
[Bhatia and Thornton (1970, 1971)]. 

The expression for the scattered intensity 
for X-rays or neutrons now involves, in general, all 
three structure factors (all RCH or NGO etc.) 4 
and experimental determination of them involves 
varying the coefficients in the expression for 


scattering. Where practical, this can be done for 
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neutron scattering by using different isotopes of one 
of the components in the mixture. The actual deter- 
mination of the structure factors has been achieved 
only for a few cases [Enderby, North and Egelstaff 
(1966), Page and Mika (1971)]. The theoretical cal- 
culations of ela) for a binary mixture are mainly 
those using the hard sphere model in the Percus-Yevick 
approximation [Lebowitz (1964), Ashcroft and Langreth 
(1967), Enderby and North (1968)]. 

The long wavelength (q>0) limits of the 
structure factors are naturally hard to determine 
from scattering experiments, and one must turn to 
other considerations to obtain information on them. 
Recently Bhatia and Thornton (1970) have shown that 
the Ay ane of the structure factors may be 
determined from a knowledge of three thermodynamic 
quantities: the compressibility, the volume as a 
function of composition and Sac (0) (the mean square 
fluctuation in concentration). The latter quantity 
(Sag (0)) depends on the second derivative of the free 
energy of mixing with respect to the concentration 
and, as we shall see later, reflects) rather sensitively 
the nature of the interatomic interactions in the 
mixture. 

In Chapter 3 of the thesis, we investigate 


the behaviour of Sac (9) and other structure factors 
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(in q>0O limit), using the well known theoretical model 
of regular solutions (in the zeroth approximation) . 
This model is applicable to systems in which (a) the 
two types of atoms are of roughly the same size, and 
(b) the differences in pairwise interactions between 
the atoms is small. These conditions are approximately 
fulfilled in many systems, e.g. Na-K. Detailed calcu- 
lations are given for the Na-K system and are found to 
be in good agreement with experiment [McAlister and 
Turner (1972)]. This chapter also discusses qualita- 


tively the behaviour of S_.(0) when the two types of 


CC 
atoms have significantly different size. 

There are, however, a number of systems (e.g. 
K-Hg, Mg-Bi, Ag-Al, Tl1-Te, etc.) where the differences 
in the interactions between the two species is not 
small. This is evidenced by the fact, amongst others 
_ (magnitude of molar free energy of mixing relative to 
RT, ©tc.)=-" that an the solid’ phase they form inter- 
metallic compounds at one or more well-defined chemical 
compositions. The main body of the thesis (Chapters 4 
and 5) is devoted to a discussion of the concentration 
dependence of Sac (0) and other thermodynamic properties 
of such compound-forming systems. 

In Chapter 4 we first describe a theoretical 


model to obtain an expression for the free energy G 


of the compound-forming molten systems. The model 
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essentially consists of assuming the existence of appro- 
priate chemical complexes Chak in chemical equilibrium 
with one another; here u,v are small non-zero integers. 
Although a general formulation is given, the actual 
calculations can be made only under certain simplifying 
approximations. In the remainder of the chapter we 
examine, in a simple approximation, the general features 
of Sac (9) and of partial structure factors for different 
values of uw and v. 

Next (Chapter 5), we discuss two higher appro- 
ximations to the free energy and calculate therefrom 
Sac (0) and other thermodynamic quantities. It is shown 
by taking several examples that the model gives quanti- 
tative agreement with experiment. 

Finally, in Chapter 6, an approximate model is 
used to discuss some qualitative features of the struc- 


ture factors of weakly interacting systems, away from 


the long wavelength limit (i.e. q>0). 


CHAPTER 2 


STRUCTURE FACTORS AND THEIR LONG WAVELENGTH LIMITS 


As mentioned earlier, there are several equiva- 
lent sets of structure factors for binary mixtures used 
in the literature. In this chapter we define three sets 
and give, for them, certain results which will be found 


useful in later chapters. 


i201 Definitions of*the Structure Factors 


The most commonly used set of structure factors 
are perhaps the partial interference functions aya (a) of 


Faber and Ziman (1965) defined by (a,8 = 1,2) 
tebe aie ig: Pg ()- Dare (22%) 


where Tyg lt) = Gp, (r) is the pair distribution function 

defined as the probability, normalized to unity for 

fatgerr, Of finding an atom of type 6-in a unit volume 

at a radius r from the centre of a type o atom. 

Because of the spherical symmetry of Fug (rv) in an 

#SOtropic- Liquid, “the Ay g (a) depends, of course, only 

on the magnitude of q (and not on its direction). 
Another set of structure factors Sug (@ 

which are frequently used, are defined as [Ashcroft 


and Langreth (1967), Enderby and North (1968) ] 
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where Cn is the concentration of a-type atoms in the 


mixture; TE Ny and N, denote the numbers of the two 


types of atoms respectively, N = Nj+ N, and? 
N N 
lic ee E 
Comparing with (2.1) and (2.2) one has 
al 
= oe: , 
Sy (a) = one + (c\cg) Lay, (a) Eliya (2.4) 


Finally, a set which we find particularly useful 
are the number-concentration structure factors Syn (GQ) + 
Sac (a) and Suc (@) [Bhatia and Thornton (1970)]. These 
are defined in the following manner. 

Let swe) (a = 1,2) denote the local number 
density of a-type atoms and let their mean density be 
denoted by ny (=N,/V) « Then the deviation én (z) from 
the mean is: 

én (r)=n (r)-n.=-n. + = Sir — Ro) ey be 5) 

a Q a 0, j=1 4 
where RO is the position: of the he atom of type a. 


Making the Fourier expansion 
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én (x) = = l N (ae ee (2.6) 


+ tae. 3 
Ny (@) = Je én, (x) d z 
No ig-R? 
= e ee oe (2.7) 
j=1 qe 


The reality of én, (r) requires that Ni (qa) = Ny (-q). 
Pah ; > 
Similarly, if N(q) denotes the Fourier transform of 


the local deviation én (r) in the total number density 


n(z) = n, (fF) + sy Gel then 
‘: ig-R% 
N(q) = e J _~ NG : 
Jeu qr 
= N,(q) + NJ(q) . LO, 


We next define the local deviation from the 


mean concentration c (=c,) by 


éc(r) = [ (1-c) én, (£) ~cén, (r)] : (2.9) 


Z2I1< 


so othat Lf én, and én, change in proportion to their 


respective mean concentrations, namely, c and (l-c), 


then $c (r) =20, assit should. ff wermakesthe Pourrer 
expansion 
> > 
6c (x) = J c(qye 9" (2210) 
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The number-concentration structure factors are. 


now defined by 


Se (a) = x <N* (q) N(q) > (22 12a) 
Sac (Gq) = N<c* (q)C (q) > (2.12b) 
Suc (GQ) = 5 <N*(G)C(G) + NIG) C*(G@)> (2.12c) 


where <...> represents thermal (ensemble) average. Like 
By (a) and Sg (ds Syn 6D) Gite.) for an 2setropic: t1dquiadG 
are also independent of the direction of q. 

To see how these different sets of structure 
factors are related to one another, one has only to 
note that the pair distribution function Tug (v) LS, by 
definition, 

> > > > 
+r) - > 
4 : <n (r,) [ng (ry .) 6g 0(r)] 
(r) = 7 d Wl. ial ite PSO ie 


g 
aB = —- ce 
<n, (r,) ><n, (rotr) > 


(2.13) 
Using for n(Z), etc. the expansion (2.5) and noting 
that <n, (r)> = nye it may be seen that Sg (a) may also 
be written as 


* 


1 a > > > 
S49 (a) Sf ee re <N* (GQ) N, (gq) +N, (a) Ne (a) > : (2.14) 
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Sincertrom: (2). 8 )ander(2% 11) 


N, (4) = cN(q) +NC(q) , N,(q) = (1-c)N(G) - NC(Q) (2215) 


one can immediately write down the S 4g (a) in terms of 
Sx (2) etc., and hence using (2.4) also the Aug (a) in 
terms of Sx (2) etc. We give here explicitly the 
relations between a yg (a) and Sx (2) etc. [Bhatia and 


Thornton (1970) ] 


2 5 ee. 
ve lana is (q) a C1 Syn (q) + 2C)5yc (q) Le Sac (q) one eS: (2.16a) 
Bee ace tiem = 2c. Sg) +S (qhecre (2.16b) 

979904 2°nn $2 2°nc ‘4 cc 1-2 : 


C4 O94 9 (4) = Cy COS (a) + (eQ~e)) Sy (a) — Sag la) + cy C2 


(2.16c) 
and conversely 
BU ay clas KG). b con, Ala) lc. ceva (a) (aia) 
OS AIR os bo a te SP wat © a Ieeeryses Lee : 
Syic (q) = C1Co [c,) (ayy (q) “a0 (q)) “Co (a5 (q) “a1> ic )uelneh.2 sega 


Sac (a) = cc, [1+ C1 Cy (az, (a) + 55 (a) - 2a,5(q))] see Qe te) 
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2.2 Sum Rules and Constraints 


It is next of interest to give the sum rules 
satisfied by the various structure factors defined 
above. 

First consider the sum rules for the number- 
concentration structure factors [Bhatia and Thornton 


CLOZO) hi ofROr Syn () we have from (2.12a) and (2.8) 
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eh) vo: Sat 2 
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1: iq r >> + 3 
= 1+5 |e ""< } 6(r-R_+R_)>d-r - N6 ; 
N mén We pearl q;o 
Hence 
: Nee | Or med Se ste (2.18) 
yas (q)-11=- B+% [ott < J oR ieee 
V. ) NN Vie oN mZzn my 1 
where we have used the relation 
> > 
i y Ra Gee L 5 pet alan Ss Soa) } (2.19) 
Me ee (27) 


The integral on the right-hand side of (2.18) is zero 

since R_ and R_ are positions of two distinct particles. 
m n 

Hence 
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wae ty a V 


' 4 
i 
Ay 
o 
{ 
\ i i 
{ ; fae 
a fj 
7 
a | y 
‘ j 
’ 
' 
i % 
} 
' " 
bAt \ i 
; k 
1 
y 
r 
x 
+ 
int ‘ a 
bi j } fi 
f 
4 é 
y i e 
aie wl 
4, OG hi 
i 
a 
: / 
i 
Ce why j 
ve aa 
Al el i 
My t \ 4 ¢ Me, Viki 
» tte A Se 
i] i@ y 
of F 


“A 


pl 


q Ay! as 


aati af sda 


“ee i: eu ves ies J 


ia 
i 
hy 
{ 
} + 
4h 
’ i 
{ 
: ‘ 
0 i 
‘ 
BT ‘ 1 
4 J 
ye f ‘ 
Ps 
Q 
in! LP a 
\! 
‘ 
ut 4 : 
r 
f ; 
j 
= 
A 
aaled Urata 4 
rie hl ee VA 
4 
Wi uf 
ih 
ie ¢ i ne 
i yt iy '' nH 
i eh 
Fy W ¢ 
Rat Cae ee 
A ey ee rn tat 


es 


r} i a: meer 


- oo MY het 


tod a (eM wig 


i 
ured 


12 


or, in the integral form 
Z aa 3) N 
ts | tsyy =e) dg = <= yo (222 a) 
(27) 
Similarly one may show that 


ee (q) - c(1-c) ]a°q = 0 (2.21b) 


[yo (a°q =o. (2.21c) 


The corresponding sum rules for a yg (a) were 
given by Enderby et. al. (1966). They may be immediately 


obtained tiromie(24:21)9 and: ((2.16;) ase«6o,8) =i 4,2) 
atl ta (Gq) wl deg =.= = Ones 
3 Fae SF Tires WAAR : 


for all the partial interference functions yg la)- 
Using (2.4), the corresponding sum rules for Sag 6d) 
are obvious and will be omitted here. 

We observe from (2.2la) and (2.22) that SiG) 
satisfies the same sum rule as avg (a) - in=Concerase, 
Phe seumeriesi (7.21) sands (2.210). tor Sea) and 
Ba) respectively, are different from that of oc) 
and Aug (d) « These sum rules clearly imply that, for 
large q, all three agg (4) and Soy (4) oscillate about 
unity, whereas Soc (4) oscillates about c(l-c) and Sue (A) 


about zero. 


Rae 


ibis “ae 


~~ 
LT 
uy'ye 
ip } 
Ws 
' ; 
te / ; Re, on 
iT 
le) 
Phi : i 1 ‘ 
{ 
<- f 
i | 
ie 
are wo 
n oe 
n a pe: ce it 
: fie . o 
‘ m 218 ; Vir ian 
’ nt ‘i ay 
| pn 
ev} 
; 1 ty | ve 
ri iT 
Ht 
a 
“A 
t I i \ ‘ } . 
7 | A 
fe J j =f 
i hit 
' a | 
{ F h 
+. } * 
: aoe = | i 
) 
ie a, 
/ 
Ae 
} j 
= us Cie Cee 


, iy of ite ai Aare a ; 4 anh ae 
y d et wa Le ' § fet iy | i & aucenl wie by ON i “ ei pid a ve af 
1 Bid Ga \ J hay 
} : : mau 


is , hake 


. ie a0 1 ‘ ut 1 _ a ul voy / 
‘'t r ; S r en T - * ee 2, eet) 
, a & j ‘ ad LO, ited am ha vie LJ “hos »} f il ‘ay 3 


Ay ah a ee | 
i ; ae \ 
i ; 
) f , j : nara ; 
) is eet f 7? thea we ony 
Co 
; i ; we y aa 
ot / 3 ray ae 
¢ j j ; inti ‘4 j ‘7 my ‘a S , sy 


€ Aik 


Apart from the sum rules, each set of structure 
factors satisfies certain additional constraints. For 
the number-concentration structure factors, two of these 


follow from the definitions, namely that 
Syn (a) 2 9 and Soca) 20. (2523) 


A further constraint can be deduced from the expression 
for the scattered intensity (in the weak scattering 
approximation)’, Lt W, (q) and We (q) are the atom form 
factors (for electron, neutron or X-ray scattering), 


then the scattered intensity is essentially given by 


1() = 18(G) 17s, () + LaWw(G) 17S, (4) +2 L(G) ] [AW(G) 1846 (@) 


(25524) 


where 


W(q) = cw, (q) + (1-c) W, (q) and AW (q) = W, (a) - We (q) é 


Coane. 


If this expression is to be non-negative for all W, (q) 
and W.(q), then in addition. to Sun (@) > 0 and S.o(q) 29, 


we must have 
BOOT Peres (ay eo (2.26) 
nw ‘2? "cc NC as 


=o A . 
in tezmms of the agg (a) » the expression for I(q) is given 


by [Faber and Ziman (1965) ] 
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I (q) = ¢W} (q) + c,W5 (a) + ctw (q) [a,, (@)-1) + 


Z2N2 
+ anon ila, Sayles 2c, CW, (q)W, (q) [a, 5 (q)-1] ’ 
C227) 


from which one may deduce that the Ang (a) must satisfy 


the constraints [Enderby et. al. (1966), Faber (1972)], 


c,4,, (a) + Cy 2 0 
Coa, (q) +c) 2 0 (2.28) 
(ca, (a) #e5) (cya59 (q)+e,) 2 ¢,¢5La,5(q)-1]" . 


We note that individually the yg (a) may be either posi- 


tive or negative. 


2.3 Long Wavelength Limits of the Structure Factors 


The properties of the structure factors dis- 
cussed in the last section, namely the sum rules and the 
asymptotic behaviour for large q, are independent of the 
individual characteristics’ of'a system: ‘in contrast, 
the long wavelength (q*0) limits of the structure 
factors depend on the thermodynamic properties of the 
mixture. This thermodynamic connection is perhaps most 
transparent for the number-concentration structure 
factors and we consider them first. 

It will be recalled that Syn (@) describes 


the correlation between the particle density and 
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particle density, irrespective of the concentrations. 
For a fluid, it then represents, in the qo70) Limit, 
simply the mean square fluctuations in the total 


number of particles in the volume V of the medium 


2 
S OLN ie 
Syn (4 > 9) = Syn 69) ey gr : (2229) 
Similarly 
x 2 ra 
Sac (0) = iN< (AG): Se, Suc (9) = <ANAc> (2°53 0)) 


where < (Ac) 7> is the mean square fluctuation in the con- 


centration, Ac being defined by 


dpe! a. 
Ac = N [(1-c) AN, cAN, ] 


and <ANAc> is the correlation between the two fluctuations 
Ac and AN = AN, + AN, - 

The expressions for the Syn 69) etc. may then be 
obtained straightforwardly from the fluctuation theory 


and one finds [Bhatia and Thornton (1970)], 


ban dc B 
Sac (9) = Sac (a> 9) = on (2331) 
=| 
T,P,N 
Z 
- Lhe = 232 
Syy (0) = Syy (a> 0) = 8 + 8°Soc (0) (2.32) 
_ = aad PA SE 
Suc(0) = Syql4>0) = -88¢¢ (0) 
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Here G is the Gibbs free energy, P is the pressure and 


6 is the dilatation factor 


V V +- V Ui e 


where va and V5 are the partial molar volumes of the 


two species in the mixture. Finally, 


Pia 
] — V k RT kp 25-35.) 
where Kp = 7 7(3 is the isothermal compressibi- 
V\oP)T,c,N 


lity at constant composition. 

If the partial molar volumes of the two species 
are the same, then §6 = O and the fluctuations in number 
density are independent of those in concentration 
(Siig (0) = 0), as might be expected intuitively. Further, 
Sx (9) Ts'then O, like the expression (1.1) for the 
Seructtre factor Of a pure liquid: 

The expressions for the long wavelength limit 
of ene other structure factors ee Pre eel follow 
immediately from the relations between these and Syn (GQ) + 
etc. given in §2:1. Substituting (2.31)=(2.33) into 


(2.16), we have (c) =c, Co = (l-c) ) 


Cc 
2 ae 2 
BOM one oh =o 4 (ee 8) 854.0) 
amas Cy Cc) Ce 
Cc 
is i: 2 
2 By 
(ee ese +S) S(O) (2.36) 
ayo (0) = 6 + 1- an = cc 


16 


Finally using (2.36) into (2.4), we have 


2 


= 1 
S10) = Se eSaSac (0) Gero) 


2 


vie: uE 
Gee yy yaa co gic, © 


eee o) = (61 ©) "8 ~ (e429) “Soq (0) (Z-- 6) (= +8)" (28 3a) 


We thus see that the concentration dependence 
of the q>0 limit of the various structure factors are 
known if we know the volume V, the compressibility Km 


and the concentration fluctuations Sac (9) as functions 


of composition. The volume and k though not easily 


Tv? 
amenable to calculation, are directly measurable quan- 
Pit Lless...For Sac (9) we require the Gibbs free energy G 


as a function of c. We recall that G for a binary 


mixture can be quite generally written in the form 


eee N tcG{? (P, 1) + (1-c) G)°) (P, 7) + Gy (2.38) 


10 g f°) and 


a 


ae are the molar Gibbs free energies of the two pure 


where W= N/N, (N. is Avogadro's number 


ST nnn nn ad UE EE ESEEnEEIeEE 


: In conformity with the standard thermodynamic usage 
we shall henceforth express the extensivity of G in 


terms of gm moles (MW) of the mixture rather than the 


total number of atoms (N). 
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Species (at the pressure and temperature of the mixture) 


and Gu is the free energy of mixing (arising from the 


entropy and energy of mixing). Since Ger 


and asc) are 
by definition independent of concentration, we have 


Promecc. 3/7) and (2.31) 


NRT 


2 
3 Gur 


0) (2.39) 


cc | 


ac 
dir Sy 
where R = NikKp is the gas constant. Thus we have to know 
Cu as a function of the composition. 


Theoretically, expressions for Gu are available 
in the literature for several models of the mixtures, 
e.g. the ideal solution, the regular solutions in zeroth 
and quasi-chemical approximations, athermal mixtures, etc. 
For reviews see, for example, Guggenheim (1952), Rowlinson 
(1969), Prigogine (1957). We examine the behaviour of 
the structure factors on the basis of some of these 
models in the next chapter. As already mentioned, the 
existing models do not in general suffice for systems 
which have, in the solid phase, a strong tendency to form 
compounds. A calculation of Gu and hence of Sac (0) etc. 
for the compound forming systems is described in Chapters 


4 and 5. 


Sac (0) can be calculated from the experimental 
data also, since the concentration dependence of Gy, is 
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tabulated [for example, in Hultgren et. al. (1963) ] 
for many molten alloys. Experimental Gu is inferred 
from the heat of mixing data as a function of tempera- 
ture or, what is more common, from the activity data. 


The activity ay of component 1 in the mixture is defined 


by 
Rpetn sake pees - clo) 
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dG 
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where A Nis Weis the number of gm moles of component 1 
in the mixture, and N’5 = W(l-c). The activity a, of 
component 2 is defined similarly. When the activities 
aceuknown, it is actually not necessary to first inte- 
grate them to obtain Gur since Sac (0) is related to ay 
Or a, by [McAlister and Turner (1972), Ichikawa and 


Thompson (1973) ] 


as (1-c) so, 2 Raa gil te (2 2Ans} 
Soc (0) ~“d &n ay d &n a5 , 
fe Ko. dc 


as may be readily verified from (2.39 rand. (2.4.00, 
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CHAPTER 3 


CONCENTRATION FLUCTUATIONS AND PARTIAL STRUCTURE 


FACTORS IN IDEAL AND WEAKLY INTERACTING MIXTURES 


So ineroduction 


The thermodynamic relations given in the pre- 
ceding chapter make it possible to Carry out a system- 
atic examination of the behaviour of the long wavelength 
limits of the structure factors for various types of 
mixtures. The present chapter is concerned with the 


study of concentration fluctuations S and partial 


CC 
structure Fy eal ais for ideal (non-interacting) and 
weakly interacting mixtures. Prior to about 1970, it 

was generally heid that the so-called substitutional 
hypothesis for a,,(q) might nee be an unreasonable 
approximation. In this approximation ane are regarded 
as nearly independent of concentration and A, 1 (A) #agg (a) = 
Ajo (qd). In the paper giving the aforementioned 

pee 29 tO Wank alley. Final $e 3 Aes oS SC Ars 


t In this chapter and Chapters 4 and 5, we shall be 


concerned with only the long wavelength limits of the 


structure factors. Hence we make the conventions: the 
words "in the long wavelength limit" are to be under- 
stood unless stated otherwise; also, the argument (0) 


in Sx 69) is deleted so that Sun= Sn 69) # Se ea etc. 
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thermodynamic relations for the structure factors, 
Bhatia and Thornton (1970) showed, by an example, 
that while the substitutional hypothesis may well 
PemreoacOnable for high q, at is far) from justified 
in the q>0 limits. Subsequently McAlister and Turner 
(1972), using thermodynamic data, obtained are fOr 
Na-K and K-Hg alloys and found them to be highly con- 
centration dependent. It was therefore felt desirable 
to examine the behaviour of 5 on the basis of some 
known theoretical models for mixtures. 

The plan of the chapter is as follows." In 
S552 , Sac and oo are described for the ideal solutions. 
In 83.3 we derive, for later use, some general formulae 
£Or as arientne- low concentration Lime. °n oo. 4, ene 
regular solution model (in zeroth approximation - i.e. 
for weak interaction) is described and the expressions 
Gk Sac and other structure. factors are derived. 
Section 3.5 is devoted to the application of these 
formulae to Na-K alloy. Finally in 83.6, some remarks 
are made on concentration fluctuations in mixtures where 
the difference between the sizes of the atoms is so 
large Prat they cannot be regarded as regular solutions. 
eee i ae Pee eee 


i The work on partial structure factors presented here 


was published in Bhatia, Hargrove anduMaren (L973) 
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3.2 Ideal Solution 
Sie fek Gyr Km, 6 and Scc 


The simplest type of mixture is the ideal solu- 


tion for which the free energy of mixing G,, is just (-T) 


M 
times the entropy of random mixing, so that 


(Gy) 5 geal = WRT [c Rnc+(l-c) & (l-c) ] Cs.) 


and hence from (2.38) 


erage ey! (lacie oie venle tn.c+ ll-cien (lecy i. 


1 2 
(332) 
From (3.2) we have for cre volume V of the mixture 
(oO) (o) 
dG dG 
dG at 2 
V= [3s -w(c(—$5) + (1l-c) (—5—) 
dP T,N,c dP . dP T 
= wtov{) + (1-c)v§$°) eee 


where sa and ve are the molar volumes of the two 
pure species (at the same temperature and pressure as 
that of thesmixture). Eqn. (3.3) “expresses the well 
known Vegard's law regarding the linear variation of 


volume with concentration. Similarly 


2 
‘3 OV Se ean | oe (3.4) 

Veg = - [Fr (25) 

VN, Tc oP 1 Ag MAT. 
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where 
(o) 
oft = - [ 1 emt ; 
fh y (0) oP 7 


is the isothermal compressibility of the pure species 
i. We note that VK and not Km varies linearly with 
efor the ideai..solution:; 

Using (3.3) "im (20:34), the expression- for 


the dilatation factor 6 becomes 


(o) (oO) 
ee = nt itat ata 
ies + (2 a (o) (0) : (3.6) 
(ey JOR - 


tBeavaLlation, Of oO *wethlc according to) (3.6), 1s) i1Ulus= 
EGatedsnef1q..3.1 for two) values. of Vogue 
Finally substituting (3.1) in’ (2239) "we have 


for the concentration fluctuations 


as = a7 
Sac = c(l-‘c) , ( ) 


where the superscript 'id' has been added to signify 
that this is the expression for Sac for an ideal solu- 


tion. 


Bee 2 Partial Structure Pactors 


With Ss given by (3.7), the expressions (2. 36) 


CC 


fox ai; become (c) =c, Co = 1l-c) 
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Variation of dilatation factor § with 
concentration for a mixture obeying 
Vegard's law. B= vito) sy fo) =e 
ral eta yoke 1 
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Ae as 2 
BR SOLON Cec G 8s 
22 1 4 Regs) 
a = 6 + §(c,-c.) + cc §2 (3.8) 
12 1 2 i aoe) ’ ° 


we observe that if the mixture, in addition to being 
ideal, is such that 6 = 0 (i.e. the two pure species 
have the same molar volumes), then the substitutional 


hypothesis is valid, that is, 


apy" G20 73 8 (3.9) 


Usually 6 is no more than 0.03-0.05 at the melting point 
G@tfan alloy and) in Figs’ 3.2 and 3.3 weisketch, taking 
0= 0, the variation of ai, with c for the two sets of 


Segivenainrrig. 3% 1 


3.3 Partial Structure Factors in the Low Concentration 
Limit 
Before considering the behaviour of the partial 
structure factors for a non-ideal mixture on the basis 
OL a model, it is useful to consider their general form 


at small concentrations. 


In the c>0 limit, the mixture approaches the 


pure species 2. Since the leading term in Sac: for 
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small c, is always! c, it follows from expressions 


(2-36) “for aa4 that as c>0 (c=c)), 


Ole i, AG Ones 0 re) 
and 
OT ar Oss Oma " 
ail = 6 26 Beco se OA Nk G3r510)) 


where the superscript 'o' means the value of the corres- 
ponding quantity is evaluated at c = 0. The expressions 
for abs and ano are the same as obtained by McAlister 
and Turner (1972), but their expression for ooh is 

ari = 8° -26°, which is true only for an ideal solution 


for which See) +9 2)e— 00% * Similar liv,e as: cy sch 
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For an ideal solution, 


fe) aes ol em 
On SIB=iLl pi oO San Hey 2 (o) (342) 


ane, acroms (3.10) "and (5.21); 


uy Ponsa, dilute, solution: .(c.<<-1)% Gur has the form 


v 
= = eee 
Gy= WRT ic Snet a,c +a5c ] 


the leading term (c &2n c) resulting from entropy of mix- 


ing, must always be there. Hence, using (2.39) one”has 
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On tak CO ar al - 
apy = O° = (By ay = ert (a-erty s . (3213) 
Biso, from (3.8) and (3212) 
a = 6 + aoe 1) a pace 
22 or c 1 


and 


0] 
2 


peers + peace se = 15 for (l=-c)e<cl.” (3.14) 


TeuLs interesting to mote that if we set 6 = 0 "in °(3.13) 
and (3.14) and identify the volume ratio 8 as (Crveccaner 
where OF and O05 are the hard sphere diameters for the 
two types of atoms, then (3.13) and (3.14) reduce to 
those obtained by Faber (1972) using the hard sphere 
Emedel. [Note the-ditrerence in notation and also that 
there is a misprint in Eqns. (6.17) and (6.18) of Faber 


(1972)]. Finally we see from (2.36) that the combination 


aiit Ano 7 2a)5 is independent of 6 and 6 and is given 
by 
Sipnjtelec) 
CC 
a,,ta.,,- 2a,, = ——|— S (Sco) 
11 2} 2 fe ame ne 


This is zero for an ideal solution, and a non-zero 
value of it is a direct indication of non-ideality of 
the solution. If we represent, for small c, the depar- 


ture from ideality by 


2 
Sh esac Se ie eee (see) 
then 
e Seen 2a ce st : cp, 
Q41 + 4822 ~ 2472 ( ) 
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We shall see presently that for regular solutions (in 
the zeroth approximation), W is a direct measure of 


an interchange energy in suitable units. 


3.4 Regular Solutions in the Zeroth Approximation 


3.4.1 Expression for the Free Energy 


For a solution to be ideal, two conditions 
have to be satisfied: (a) the difference in the sizes 
of the two types of atoms must be small - in terms of 
the volume ratio 8 or (eso) it is usually considered 
Sufficient if > < 6 < 2 [see, for example, Guggenheim 
(1952)]; (b) the difference win the pairwise interac- 
tion energies between the atoms must be zero (see below). 
In the regular solution theory the assumption (a) is 
retained, but the difference w is considered to be non- 
zero. | 

When w # 0, only approximate expressions for 
the free energy of mixing Gy can be derived. The sim- 
plest approximation to Gui which is valid when w can 
be considered to be small (see below), is usually 
referred to as the "zeroth or Bragg-Williams approxima- 
tion" (Guggenheim, 1952). The expression for the free 
energy of a regular solution in the "zeroth approxima- 


tion" is usually derived in the following manner (e.g. 


see Guggenheim, 1952). 
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One assumes that the interaction is signifi- 
cant only between atoms which are nearest neighbours 
of one another. Let Xanr Xpp and Xap denote, respec- 
tively, the interaction energies between nearest neigh- 
bour A-A, B-B and A-B atom pairs’, and let the number 
of nearest neighbours of each atom be Z. Then if Naar 
Nop and Nap denote the number of nearest neighbour A-A, 
B-B and A-B pairs in the mixture, the energy E of the 


solution is 


B= (3518) 


NaaXaa * NgpXpp t NapXap 


RIN, aN =5NZ = 5NN2 , (3.19) 


N being again the total number of atoms in the solution. 
If we now assume that the two types of atoms are dis- 
tributed completely at random, then of the Z neighbours 
of a given atom, a fraction c will) be A atoms and) (1-c) 


B atoms. Hence (c= Cy) 


=2(i 2 (3.20) 
Nap 25 NZ je Cul=C) pec. 
iyo ee eA ee eee. Bee Fe eee 
t For convenience of writing two types of atoms will 


sometimes (as here) be also denoted by A and B with 


A=species 1, B=:species 2. 
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Substreuting =( 32 20) am (3.18). and rearranging terms 


tee a 2 2 
E = pare Xan t (l-c) Xpp t 2c (1-c) xy_] 
af 
= Wwe(l-c) + 7 (WNL Z) [cXan + (1-c) x, ,] (3. 22) 
where 


w= (5 NLZ) (2Xap - Xan — Xpp) - (3.22) 
Since the terms in the square brackets are the appro- 
priate linear (in c) combination of the characteristic 
energies of the pure constituents, the energy of mixing 


Ey is simply 


Ey, = Wac(i-c) . (3.23) 
Finally, since the atoms have been assumed to be at 


random, the contribution to Cur from the entropy of 


° e e ° . oF 
mixing is the same as for an ideal solution. Hence 


t We note that this expression £OG Gur is the same as 
that for the conformal solution node lalenateren caine 
(1951)]. Here one assumes that the properties of the 
mixture can be described by small perturbations from a 
suitably chosen standard reference liquid (at the same 
temperature and pressure as the mixture). The word 
“conformal” refers to the fact that the interatomic 
potentials for the different pairs are assumed to have 
the same functional form. The two models are valid 


essentially under the same conditions, namely the assump- 


tion (a) and smallness of w. 
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G,, = E,, + (G_) 
M M M ideal 


=WRT[c 2n c+ (l-c) $n (1-c) ] +wwe (1-c) (3524) 
and 


(o) 


G =WcG 


ie Pati) G 


+Wuwe(l-c) , (3e25)) 


which are the desired expressions for Gur and. G for a 
regular solution in the zeroth approximation. 

The quantity w, assumed to be concentration 
independent, is called the interchange energy. If we 
replace an AA pair and a BB pair by two AB pairs, 
then the energy of the mixture is increased by 
[w/ (5 N.2)1. Clearly if w is positive the formation 
of like atom (nearest neighbour) pairs is energetically 
preferred over that of unlike atom pairs. The opposite 
is the case if w is negative. This tendency of an atom 
to prefer one type of atom over another as its nearest 
neighbour was neglected in the above treatment when we 
assumed the atoms to be distributed completely at 
random. A somewhat better approximation is the so- 
called quasi-chemical approximation. Here it is found 
that the various types of pairs obey a sort of chemical 


equilibrium relation, namely, 


3 +NVRT [c ln ct (1-c) £n(1=—c) J 
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W 
po 4 
TARVER Beira tals! 4) 2u/(2RT) 
a oR e = e ? (3.26) 
25 AB 
dnecontrast to. (from .(38.20) ) 
iN. —N 
AA BB iw 
inne = 1 (3227) 
2 AB 


in the "Zeroth approximation". Since the quasi-chemical 
approximation is itself approximate and quite gene! 
we shall not pursue it here. 

A comparison of (3.26) and (3.27) shows that 
the "zeroth approximation" should be valid if 
|2w/(ZRT) |<< 1. If we say, somewhat arbitrarily, 
that the approximation should be valid for |2w/ZRT| < = 


then since Z%8, 10, we get a rough criterion for the 


eZeEvroOth approximation, tO be reasonable was 


Ww 
728s eee = (34218) 
Since (Suin given by expression (3.24), is maximum at 
c= the, condition (3. 28)Rimpliespinar 
G 
M it 
~ ~ = 5). x sar Ab) 
Pe WRI Oe on Nate Cc x) ( ) 


We shall need (3.28) and (3.29) in our discussions 


later. 
pera tian) cree rrees 3 Se 


t For reviews see, for example, Guggenheim (1952) and 


Prigogine (1957). 
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3.4.2 Expressions for V, 0, Km and Sec 


With G given by (3.25), one obtains for the 


volume, § and Ke the expressions: 


= wt tovi?) + (1-e) vf) + atv! e(1-c)] (3.30) 


6=2(2| = ie AS a a'vi°) (1-26) ay 
ae eae : 
Vide A as lar ev? + (1-e)v5° +a vt?) (1c) 


Ae te ee? 4 qty tO) Mara 


(3.32) 


where a' and a" are dimensionless quantities proportional 


to the first and second derivatives of the interchange 


energy w with respect to pressure 


a 
} i dw " as d* w 
a — ane Qa =- ey e G53) 
(o) [Se d O or| | 
vy i Ns Kn dP T 


If w is independent of pressure, V and Vk are linear 
in concentration as for an ideal solution. 
Using o(S 724) inac2)..39)one nasitor sehneycon= 


centration fluctuations: 


S = ee : (anon) 
CC : 5) hoe ee (l=c) 

18 RT 

9 Jw,T,P 
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O 0-5 
C 


3.4. Variation of Scc with concentration for 


regular solutions in the zeroth approxi- 


mation; w/RT = 1.1, —— w/RT=0.5, 
ecoeee W/RT=0 (ideal solution), --- w/RT= 
0.5, —— - —- =W/RT = -1.0, @ data for 


Na-K from Hultgren et al. (1963). 
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We observe that Sac is symmetric about c = > and 
lies above the ideal value ee = (1=c))] pasty urs 
positive-and below it af.6 +s negative. If wis 


positive, the solution has a critical point of mixing, 


the critical temperature dU nak being given by 
ee Ww 
T =5e ‘ (353 53) 


Near this point §S Cc becomes very large. Although the 


Cc 
expression (3.24) for the free energy is too simple 
to adequately aasetne partially miscible Iiquids |, and 
the thermodynamic theory of fluctuations breaks down 
near the critical point, the largeness of the concen- 
tration fluctuations near ee is of course well known 
from the phenomena of critical opalescence in light 
scattering and other scattering experiments near q = 0. 
Eig. 3.4 allustrates, the wartalionsot Sac 


with concentration for a few values of w/RT. 


324.3) Expressions) for the Partial Structure Factors 
By combining (3.34) with (2.36), one immedi- 

ately obtains 

13.15) S20 Jor ee jee) ee US SN 

iPneiasiecnesed: neare (3.28) 4 fOr |w|/RT > 2, the zeroth 


approximation is likely to be poor whether w is posi- 


tive or negative. 
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4 2 
a4, = 8 + X(6%c(1-c) - 26 (1-c) + a0 (1-c)*) 
2 
ao = 6 + X(d c(l-c) + 26c +28 e*) 
and 
Bee Sania xis ee oyeeemiene) 6s le 
12 se c(l-c) ) (3.36) 


where we have set 


1 - (Rp) c (1-c) 

If we take the limit . c;0 (c =C)), then 

Osun =» AOl: O 2W) 

aT = § 26 ap RT 

Oy 20 

a55 = Is) 
and 

fo) .e) fe) 

a1 = 6 - 6 , (34:38) 
yielding 

fo) fe) ey 2 2h) 

ail + a5 7 2aj5 = RT (3:339)) 


and showingmthat w)) in Eqns (3 3)7) ysvsimply related to 
the interchange energy #® in units of RT. The left-hand 
side is zero if w = 0, in agreement with the result in 
(So. 5), that for anyideal solution 441+ 8o27 2412 = 0. 
When w # 0, this combination of ais is given by the 


expression 
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1 - c(1-c) (24%) 


Its variation with concentration is illustrated in 


PiG.4 3.10. 


3.5 Application of Regular Solution Theory to Na-K 


3.5.1 Comparison with Experiment for Scc, V and Kp 


For the Na-K alloy the data at 300°K on Gur 


and activities are given in Hultgren et. al. (1963). 


Within experimental errors G,, is symmetric about c= 


G 
ree Peeves 
and,ac.c. = > and par 


IS) 


M 
-0.43 so that the regular 


solution theory as described above should be valid. 


The determination of S from theory just requires one 


le 


parameter, namely the interchange energy w. Choosing 
w so that Gy/V NT agrees with experimental at c = > 


one finds 


O/,RT = AL, Te = 52.9 Olen Oe ane 


With this value of w the plot of Sac against c is given 
invbig. 3.4. As described’ in $2.3, the experimental 


values of S are most readily determined by differen- 


ce 
tiating the activity data. We see from the Figure 
that the agreement between theory and experiment is 


fairly good. 
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For calculating V and Kp as a function of con- 
centration, one needs, apart from the volume and com- 
pressibility of the two pure species, the first and 
second pressure derivatives of w, i.e. a' and a" defined 
in (3.33). A theoretical calculation of these, as of w 
itself, is a difficult problem and we regard here a' and 
a" also as parameters. Using the observed volumes and 
compressibilities for pure Na and K, Figs. 3.6 and! 3.7 
illustrate the variation of V and VK with concentration 
Clror several values of ‘a* yand a™. The curve with a? = 
-0.1 for volume and the curve with a" = -0.5 for VK 
fit quite well with the respective data on these quan- 
tities due to Abowitz and Gordon (1962). Thus we see 
that for the Na-K system the departures from linearity 


are quite small. How these small departures affect 6 is 


denircted in Figy 3.8. 


S222 Partial Strucbure Factors 


Once S V and K and hence 6 and @ are 


ceil ii 
known, it is an easy matter to evaluate ais uSsings (312364).. 
From the values of V and Km for the Na-K system given 

in the preceding section, one finds that 0 [= (N/V) RTK] 
lies between 0.05 and 0.0 in the entire concentra- 
tion range and makes negligible contribution to oa et at 


most concentrations. Fig. 3.9 depicts the variation 
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----- a' =0, —-—- oa' = 0.5. 
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of ais with concentration for the parameters ety ee 


and w appropriate to Na-K system. Since as discussed 
above V and VKn vary approximately linearly with c, the 
curves-taking a' = a" = 0 are also given. We see that, 
dopexpected,, the two sets of curves lie very iclose to 
each other. 

Fig. 3.9 also gives the experimental values of 


a4 ebvained: by McAlister: and Turner (1972). Pai will 


be seen that the agreement with theory is good except 


fOr 4a Ce y nearwc C=C. 8) = -0 sand Lora (ae) 


11 ‘"?NaNa i 


near c) = it This “dirserepancy is’ due’ to the act that, 


Na 22 KK 


as discussed in 83.2, they determined aaa (value of 
aii at c= 0) .by substituting the, observed values of 


eo and 6 “in the’ formula 


oe GO a 8B (3.41a) 


: id 
which is only appropriate if Sac = San Tecate tian 
ine expression (se lo), wirch for our moder, dives [from 


(37738) 


Oa ue oO, es (3.41b) 


The values of an from msec ta )e and ?(5.41 by ™are= respec= 
tively 0.7 and=279%atel00°C.. “Theuse of sthe latter 
value for ari would remove the sudden drop near c= 0 
in the experimental curve for ayy thus bringing it 


close to the theoretical curve. Similar remarks of 
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course apply to Aao5 near c= 1, the two values of ans 


now being 


oe iearal: as 
and 
jee S ii PAA 
a50 8 + 26 AP RT GW § Pe ds ie ee 


Be2.5 Concluding Remarks 


To conclude, we see firstly that for the Na-K 
system the ans vary markedly with concentration and 
that their values are in general very different from 
ais =~ §@ expected from the substitutional hypothesis. 
Secondly, we have shown that the regular solution model 
in the zeroth approximation is capable of explaining 
allethe data,forithis system. ..Ror-K-Hg alloysalsomi(the 
other system examined by McAlister and Turner), the 
ais vary strongly with concentration. We cannot how- 


ever apply the above model to this system since the 


Magnitude of the observed free energy of mixing 
es se eat ices 0.4m, este o00sK., 


in conjunction with (3.29) shows that this is not a 
weakly interacting system. Also, the volume ratio 


B (= Tea Me = 3) is large. We shall consider this 


system in the next chapter. 
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3.6 Effect of Size Difference on Concentraction 


Fluctuations in Weakly Interacting Systems 


We conclude this chapter by briefly discussing, 
for the sake of completeness, the behaviour of concen- 
tration fluctuations Sac in a system for which, though 
the interchange energy w is small, the difference in 
volumes of the two species is too large for the regular 
solution model to apply. For such a system a more valid 


approximation, in a Similar vein to (3.24), is Flory's 


formula (Flory (1958) ,“Guggenheim: (1952) ], 


Gy= RTI [c Qn o+ (1l-c) ¥n (1-4) ] }+{Wwo (1-6) [c+ (1l-c) /8] } 
(3542) 


where, as before, w is a suitably defined interchange 
energy, while ¢ is the concentration of species 1 by 
volume and 8 is the volume ratio. With w independent 
of pressure, as is usually assumed when using this 


approximation, 


ere. ; yao ee 
sg a TS CNT SA peor : 
ae i (1-c)V5° We2 
From (3.42) and (2.39) we have 
S = ’ - 
CO 4 cli-e)é* (1 - —2829 1 
(B-1)~ RT 


where 
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We see from (3.44) that, unlike the case of 


regular solutions (6 small), Sac (0) now also depends 


on 6. In particular, even for an athermal mixture, i.e. 


when the interchange energy w= 0, Sac is not just equal 


to c(l-c) as for a truly random (ideal) solution. The 
variation of Sac for several values of w/RT, for the 
case 8B = =) is) viiustratedmin: hig. -5..0. 

We should finally mention again that (3.42), 
like (3.24) for the regular solutions, is valid only 
for weakly interacting systems with restrictions on 
Gu Simitar to 3% 29)". Also, we. should note that 
although Flory's expression does indeed reduce, for 
B§= 11 to, the "regular. solution” expression (3.24), at 
should probably be viewed to be valid only for about 
> 1 
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Fig.ws..10 Variation: of4Scc wth? concentration 
in Flory's approximation (eq. (3.44)) 
for the volume ratio 8 = 1/3 and 
several values of w/RT (Fig. from 
Bhatia and March (1975)). 


ie | 
ay | , 
\ io y is 
f 
Minis a —_ wa ‘om mae mel waip err 
f : i ; J 
i 
f 
i tf 
{ ‘ ’ 
ie bs] 
i] i 7 
‘ " a 
t F 
‘ ' " 
Vie 
i 
5 1% : 
‘ r ’ 
a nt! r 
‘ ge atins od f 
, ya. 
1 ¥ ‘ 
f i 
‘ 
, 
4, 
i 
sae Drs 
i Bay 
A 
i \, 
% 
™" Mi 
y 1 ' 
4 . 
” ‘ 
e hel ‘ 
i } ' 
is 2 A 
5 
Sa 
- % © 
é As 
Py 
q ‘~ +, ~‘ 
ie ¥ 
ilk) ' 
eg 7 i 
"Paiute a 
il aa, ‘ 
v 
' f : () { + 
} i } eM i 
i f - ' | | ; ‘ 
} j i {bey re i ?, lee : al r % a f 
” wen + ele ties + epha patg gta wyerngat Ah ret orm at did yt ddvy me ay le mins nr anit 
ai ; Y 
i f° K wi 4 
Tqyl 
ere aii 
by | , 
7 i ‘ A 
* 5 ull Ao 
. T ‘ a | Pe is 
ce 4 . y uy ee | bas 
’ a u 1, : 7 
re hes sit } Na 
; hut y y ha i rae fae ie 
i ‘ ’ ' re aes) ie oTY ey rf M4 7 ra * id 
‘ ; Spun! | ah nes iy en a Nag: 
we a y Kad . 
‘ ‘ 7 a toh) ry? Ld eo rate iby a e ane ye A uM 
7 ‘ i ek i ; 4d) mi 
atlas SNL ie Ai Sika Dh AS), ob} Ki; 


wea, sea? os, ia ety ae: ie) a." 
te tx ; ; bial Sag aE sy 


- AL. LD i 
a ; ni y Wey : wy ) ti 7 
; 7 i a ; f “a 
Seppe: Vas ees: ae 
Ge i _ ' 4 ah ee 


; ; 7 7 7 ‘ 4 ; F on, & an iF 

> si Se oy ! (bos a 3 | aa 
‘ : yd a 

- 5 i 7 _ i 
~ AO ae 
ee | 


CHAPTER 4 


CONCENTRATION FLUCTUATIONS AND PARTIAL STRUCTURE 


FACTORS OF COMPOUND FORMING BINARY MOLTEN ALLOYS 


4.1 Introduction and Model 


In the preceding chapter, we discussed the 
behaviour of the concentration fluctuations Sac and 
other structure factors on the basis of some simple 
theoretical models of mixtures. Basically this 
discussion was applicable to systems in which the 
difference Ww in the pairwise interaction energies 
between the two species can be considered to be small 
‘(see Eqns. (3.28) and (3.29)). There is, however, an 
interesting class of molten systems - examples being 
K-Hg (already mentioned in Chapter 3), Mg-Bi, Tl-Te, 
etc. - where this difference is not small. This is 
evidenced by the fact, amongst others (magnitude of 
Gy/YV RT etc.), that in the solid phase they form inter- 
metallic compounds at one or more well defined stoichio- 
metric compositions. In this and the next chapter we 
discuss, using a phenomenological model, some charac- 
teristic features of S and other structure factors in 


CC 


; “ 
such "compound forming" solutions . 


cetperereaptchcceecenetanctteren dtm eA ET ETL TTT 


t The results described in Chapters 4 and 5 are pu- 
blished in Bhatia and Hargrove (1973, 1974), Bhatia, 


Hargrove and Thornton (1974). 
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The concentration fluctuations Sac have been 
recently evaluated from thermodynamic data for a number 
of compound forming molten alloys by several authors, 
notably, Ichikawa and Thompson (1973), Thompson (1974), 
March et. al. (1973), McAlister et, al. (1973) and 
Turner (1973) ,~and Figs. 4.2-4.4 and-5.3(a-d) exhibit 
the variation of Sac with composition for a few systems. 
A remarkable feature of these results is that, even for 
systems in which the difference in volumes of the two 
pure species is small (like Mg-Bi, T1-Te, Ag-Al), Sac 
is far from symmetric about eek, This) is¥quitesrn 
contrast to the expectation from the regular solution 
theory which in both and the zeroth and higher approxima- 
tions predicts that G,, and hence Sac be symmetric about 


M 


c= (with the usual assumption that w is independent 
Cc 


To interpret the behaviour of Sac and other 
structure factors in compound forming molten alloys, 
one has thus first to devise an appropriate expression 
fOr Gus Rifirst principle calculation (fromestatis= 
tical mechanics) of Gu is difficult for these systems 


and does not seem to be available. We therefore here 


use a phenomenological model. 
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It is assumed that the mixture consists of the 
A and B atoms and chemical associations of the type 
ALB. (u,v small non-zero integers) in chemical equili- 
brim. The requirement that, at constant temperature 
and pressure, G be a minimum then gives the equilibrium 
numbers of different chemical species in terms of their 
chemical potentials. Although a general formulation 
is given, the actual calculations in this chapter are 
carried out by making the simplifying assumption that 
the mixture of different chemical species can be con- 
Sidered as an ideal solution. For liquid (as opposed 
te gaseous) mixtures this assumption is admittedly 
approximate. However, it does make the problem consi- 
derably simpler and has been used in the past with vary- 
ing degrees of semi-quantitative success to understand 
a variety of properties - see, for example, Bent and 
Hilderbrand (1927), Barfield and Schneider (1959) and 
Darken (1967), where other references to the literature 
may be found. In the next chapter we shall consider 
improvements upon this ideal solution assumption which 


will allow a more quantitative comparison with experi- 


ment. 


t The effects of such chemical associations on Sac has 
also been considered by McAlister and Crozier (1974) 


using similar methods. 
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The plan of the present chapter is as 
follows. 

In 84.2, the model is described and on the 
basis thereof the formal expression for Sac is derived. 
In 84.3, this is specialised by making the ideal solu- 
tion assumption and is then used to discuss Sac for 
specific cases of different u and v (§4.4). For com- 
parison, Sac aS evaluated from measured thermodynamic 
data for Mg-Bi, Ag-Al and Hg-K systems are also given 
here. Finally, in 84.5, we apply these results to 


discuss the variation of partial structure factors with 


concentration. 


4.2 Expressions for G and Scc 


A distinctive aspect of the molten mixtures 
under consideration is that in their respective solid 
states, they form compounds at one or more definite 
chemical compositions. It is therefore reasonable to 
assume that in the molten state of the mixture there 
is a strong tendency for the two types of atoms A and 
B to form chemical associations of the type AB, 
(where yu, v are small non-zero integers), the precise 
numerical values of yu and v depending on the mixture 
under consideration. 

To write down a phenomenological expression 


for Gon the above basis, consider, for definiteness, 
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1 gm mole of a mixture of A and B atoms, the molar 
fraction of A atoms being c. Let at any instant these 
exist as nj, gm moles of A atoms, n5 of B atoms and ny 


gm moles of A B,, molecules, where ao = 3,4,....m, and 


Ha Yo 
Uy and Vy are again small non-zero integers. From 
the conservation of total number of A and B atoms, we 


have 


and 
m 


NS ‘ie n,= 1- es CY a el a (4.2) 

Se G. denotes the chemical potential per gm 
mole of the species, i=1,2,....m, then the differential 
change in the free energy G = ) n,G; Ls) given by (S, 
entropy) 


m 
dG = -SdT+vdP+ }) G.én, . (4.3) 
aig ees 


We can now regard nye a= 3,4,-.eeem aS order parameters 

eee eine es Ga values are determined by the require- 
ment that for equilibrium at constant T, P andc, G 

is a minimum. Hence eliminating ény and én, from (4.1) 


i = 0, we obtain 
and setting SDE Sa ’ 


HS + V S27 Gyre Ps CQ Sed (Aire ollie. (4.4) 
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EE Gir ESL, 2y.. aM, are known, Eqns. (4.4) together 
with (4.1) constitute m equations to determine the 


equilibrium values of Nyr Nor cece na at a given sf, 


Peandec. 
It is convenient to write the chemical poten- 
tial sqinetheforms (i= s),,2 ».32.6m) 
c.= 6! 4 RP en(n,/n) + RT en (4.5) 
i i a aa? t 
where Ges is -artunction of P and T only and is, just 


ehemmotar Gibbs free energy for the areas species i (at 
the pressure and temperature of the mixture) and Weg by 
definition, is the activity coefficient of the species 
i. The ¥,'Se in general, depend both on the various 
concentrations fale and Onewe and ff. )Lliethe mixture os 
different chemical species can be considered to form an 
ideal solution, then all the Yi = 1. (This assumption, 
of course, does not imply that the molten system consi- 
dered just as a mixture of A and B atoms forms an ideal 
Solnecion.)~ Substituting (4.5) ,1ntot (494) Faone has *kor 


the equilibrium values of n; the equations 
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of ideality Cy; Spee Li Zyp eee .M),. Bans. 14.6) .become 
the familiar equations for studying chemical reactions, 
Ky being known as the reaction constants. Note that, 
to avoid cumbersome notation, we have denoted in 

(4.6) the equilibrium value of a n, also by n,- In 


the equations that follow all n. Signify equilibrium 


values. 
The Gibbs free energy may now be written as 
G= ) n,G, 
ve 
= c61°) 4 (1-c) 65°) - os NG ,+ RI } ni: pena, ,7n) + 5 niA, , 


i=l i=l 
(4.8) 


where A; = RT &n Y5° Noting that Gy are independent of 
the concentration, we may then obtain by straightforward 


differentiations and use of (4.1), (4.2), (4.6) and the 


Gibbs-Duhem type of relations for Aye 
m m oA. 
v= [36] =ov{°) + (1-c)v3°)- bn oa ns [se | (4.9) 
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where v_ = (dg /9P) , and a prime denotes differen- 
Qa a Jur ie 

tiation with respect to c. Further, ye and Vice 

are the molar volumes of the pure liquids A and B 


respectively. The expressions for § and «x, may be 


x 
obtained by differentiating (4.9) with respect to c and 
P, respectively. 

In order to apply Eqns. (4.6)-(4.10) toa 
particular system, we have to specify the y,'s era Wye 
Vat As already mentioned, the simplest approximation to 


Xe results if one assumes that the mixture of different 


Chemical species M. (M. =A, B, A B,, ) is ideal, so 
a al: Hy Va 


that vall Y; = 1. This is tantamount to assuming that 
there is no interaction between the different M;. The 
next approximation to Xe would then correspond to assum- 


ing pairwise interactions w.. between the different M, 


ij 
which are sufficiently weak for the theory of regular 
solutions (in zeroth approximation) or of conformal 


Solutions, to apply. Foruthis case y; are given by 


[Longuet-Higgins (1951) ] 
Noa exp [A,/RT] 


m 
= exp{I } 


, (94/9) 945 - Te ies (ny /N) wy, 51/RT} (427A) 


with Wis =-09 tf-i-= 7. rhe use~of ie From, (4.1 )e; 
apart from adding a number of additional parameters 


(wi5)e makes the problem of solving (4.6) difficult. 
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As mentioned earlier, we shall adopt in this chapter 
the ideal solution approximation (A; = OF ¥, =), with 
the reservation that to obtain detailed quantitative 
agreement more sophisticated expressions for ne such 
as (4.11) will have to be employed. 
There are a number of binary alloy systems 

which, in the solid state, form compounds only at one 
stoichiometric composition, e.g. Mg-Bi at Sug = 3/5, 


Ag-Al at Cag = 3/4. In such cases one has to consider 
only one pair of values for yu, v. The same is appro- 
ximately true when the compounds are formed at more than 
one stoichiometric composition, provided one of these 

is considerably more stable than the others. For 
simplicity, and also since our main purpose in this 
chapter is to discuss broad features of the concentra- 
tion Pigetadtions and sof) the partial ’structure factors 
rather than consider one specific system in detail, we 
shall present here results of calculation assuming that 


only one type of chemical associations are formed in the 


molten mixture. 


4.3 Expressions for Scc etc. for Single Pair of (uey) 
oe a ee eee 


and the Ideal Solution Assumption 
It is convenient next to rewrite some of the 
above formulae for this simple case. We denote the 


single type of chemical association by ao and its 


. 


ja: tc oe 


mi ¥ ral « a 4 


eterno het Paget, piso 9 hater tess Spine) 


| ee ey HU pa. No ¥ 
ae -o f ity iy se bet Ghat aay ie i o satin h\ 
' rey ye 
i ' _ MBs j 
7 ae | Beh | Ore 
, { a" 
‘cys f HN) { ‘al esi. v7 Cee al in) 
ari , i Oh ide 
li 7 f » 
, i Mae 
\ ‘ aa ' Ww ape en 
j ‘4 i! de 5 pet lh ‘ da o L t (ree: ] 
an m Pi 
{ ae 
ARI -) 
ih ine " et it ie a 
rrig 3 mY ah 43) 6) Os ys ct} ue esi 
hey Ge i NAT mes 
; eae ha , ii MRED \ 
iy hee bb ne rT) oh I es a et ot Mc Tee AG ft ee ee 
au shy 1 te hae A ie ay eh HANS A an a ae 
ni . " pate ty at ed om” thay f 
fl re iy; r 4 i] } 
( path 7% itis 
> ro. eae t 
LEYS int) sae 
* 7 as, 
| Aiea 
A ' 1 
er | oy 
erin es i ee eee, 


' Al si eee) ae foe | i l : i hy “a ne 
t ” n Pe, ip i . i‘ ae 


is ak ve ibs 
a al soa 
Nii 1 me 


‘han hi 


; | i ! : ro 
. y bie i ae 


i ' nN } ee ACh Ce ve Te 
Ae. dg Os: i y ooh Relea) My yt Diy | Bot a Fy ive. aT, en ee 


one AA Rl Bin. mL MOR, 


i en 
rl ‘rh F 
a 2 
P aagee aie | + ey 5 ae ee 
ru! vigt Da ey a ‘ 
> aes a ia Wt 
: ne nh, id 
| , ey nee way aD i ¥ ; 
7 " ’ ie ‘ ma) Beit: i 
Ca Gee | } 
; / 4 M1 a a Wp . * 
; ' ‘i y : ’ } tl . cre bal A wit, | a wu 
rh I : i my, 
i tay ‘a ‘ iH) 
Pint ae ‘ 
; yD eran 
hall é 4 . ' é i ‘ * fs see : nul ae 
, a : ¥ I 7 
s ‘; a an ne 
vig ay ™ 


te ited 


‘ne 4 
eCg 


vit a » TUE SCS PAs te f) sie : sue yee 
h Li = ‘ ie j ie 
f ' ‘ mi rae | i j ii / 7 Thee | ie .S A 


- . tli mara a “ r 
‘ he waa ‘ a ‘ ay) 
: ’ * el q in ‘ 
< roe 
tag 5 


molar amount by n 


co 


ny=C- Wn, , N= LiraCis vn3 n=1l- (utv-1)n3. 


3° Eqns. (4.1) and (4.2) then reduce 


(4.12) 


Purther, (with Voce, "oem lo. iy Bans: (455 )= cour (4 a0) 


yield 
(n,/n)* (n5/n) ¥ 
eee ee = exp[-g/RT] =) K ’ 
WiLth 
g =-RT fn K= ven a se = on . 
(0) (0) : 
as cg, + (=a) ¢5° SPDs IY STR Rn eu an(ni/n) , 
v= ovis (1-c) vy0- nav fi Vis sRT Go) Yn K/ OP) 9 ; 
(oO) (o) ' 
6-3 (2) ON a ota 
SVvroc a (0) (Coe : 
ip ep cVy + (1-c)V, N3V 
RT / OV RT (ope mGl,) (oyna) p) 
Ae ea - Sev cid. ely dal aatre, von) 
yea ge y2t i oT Sail oP 
and " 
| S = ; (n,) eh eye 
Ge 7=L ny . 


From (43.12) 


nj= 1l- uns ; a A vn on = tT. 


(4.13) 


(4.14) 


(4.15) 


(4.16) 


(4.17) 


(4.19) 


(4.20) 


— 


ae es, ee 


Hence remembering that K is independent of concentration, 


we have from (4.13) 


: [(u + vic-plnn, 
23 = >? CO C“(4.221) 


2 2 
nyNnon, (utv-1) = Honnon, - nV) nn ,n3,-nn)n, 


Care must be taken. in using (4.21) at ec = 0 andc=H1, 
since here nN and one of ny and n5 simultaneously tend 
Lo-Zero. 

Finally in evaluating 6 one needs Soa Belcan 3 , 


wnaich is’ given by, from (42.12) and (4.13), 


-l 
on Pa sg eee ee: I heey 
IP <enp Cnet? Foe et is any, coe ete (ooze! 
Lec ah 2 3 


We observe from the above equations that the 
determination of concentration fluctuations (Sac) 
requires only the knowledge of the reaction constant K. 
An approximate value of K for a specific system may be 
inferred from (4.15) and the experimental data on the 
free energy of mixing Gu = G= eG) Slee) Gy° The 
determination of V and @ from Eqns. (4.16) and (4.18) 
requires the knowledge of v and (dv/9P) 7, respectively, 
which are related to the first and second derivatives 
of K with respect to pressure. In practice v may be 
inferred from the data on volume. 96 is usually small 
(v0.03) and may, in the absence of experimental data 


and to a first approximation, be linearly interpolated 


between its values for the two pure liquids. 
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4.4 Concentration Fluctuations Sac for Different u,v 
Ser ees ee ee eee gee ee a, Pe et cs ern 


We now use Eqns. (4.13) and (4.19) to calculate 
the concentration fluctuations Sac as a function of con- 


centration and the reaction constant K for different 


pairs of values of u and v. 


4.4.1, Limiting Cases 


If we eliminate n, andn using / (42.12). Eon. 


1 2" 
(4.13) becomes an algebraic equation of degree ut+v in 
N3- Hence, except for the case uw = 1 and v = 1, one 
has to solve (4.13) numerically and then obtain S.- 
‘Before presenting the numerical results it will be use- 


fat first to, consider icertain Jimiting :cases. 


a) K > o 


First, if K = ~, one can readily see from (4.13) 


that n, = 0 and also Ce jai, =.0.. . -Lhnen Nn) = ¢, No= l-c, 


5 
and one has from (4.19) and (4.20) 


id _ ‘ (4.23 
Sac = c(l-c) , ( ) 


where, as before, the superscript id on Sac indicates 


that the binary mixture behaves like an ideal solution 


in this case. 


DSi 50 


The solutions in which molecular formations 


occur to a marked degree form the opposite case. Here 
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i -4 


Meoeel jew tn typicalwyaiues.of 10 “)— 100. “We there= 


Poresconsider (the limiting case K = '0), although physi- 
Ca Wl yek—cany never berzeropas) it comrespondsmto anfinitely 
Pace negativerstree venergyeor mixing™ (choy Eqn. (4.15) ie 
We shall see presently that the simple expression for 
Sac for K = 0 corresponds closely to the Sac for a non- 
Zeon <:<. 1, 

With K = 90, (4.13) implies that either n, = OF 


so that yy OY It sOre ny = 0 so that nz = (k= ce 


Remembering that the maximum possible value of nN is 


a ie eeegs eee we have 
n3=C/y iebe OM CASA) ; C424") 
mm, Silce) /y fou L/ (ies ac a< el: ; (As 25) 
Pierpconcentrations c=) 0) and) c = 1 arevexcluded from the 


Vanges occ in (4424) sand o(4.25).) Here N32 = OSso4thas 
the plent ihand. side ‘of eqn... (4. 13) has aene form 0/0 4and 
it needs special consideration to which we come 


presently. 


When (4.24) is valid, one has 


De, ny = 1 Ct =tGN / 1. ; n= b= tiev=Dic7 i) (4526) 


| 
= 


Z 

' 

and (nj) /D4 = 
Hence using (4.19) and (4.20) one obtains 


ce a fu-c(ut+v) ]fu- (utv-1l)c] , O< CS TESO (420273) 


A t 
hs \ 
=) 
I 
ri 
LA 
4 
i ai i 
‘, 
{ 
, “ 
' 


i 
a 
1 om 
ti \ t } 4 
H ’ 
\ a 
\} 
Wu i 7 
Nd 
oF , by ica) 
Hatha weyers 
v * 
rn 
ve 


ibs ie sa ti 


hid 


6) 6 ag 


64 


Similarly when (4.25) is valid 


eos tic) 


ce = [v-(l-c) (utv) ] [v= (l-c) (utv-1)], —-hec<1. 


(4.28) 


Note that (4.28) may be obtained from (4.27) by replacing 
everywoere in §(422/)>c by l—-c, joby v and v byt. Since 
(4.04) was, in general, not: valid at c = 6; aisccannot 

be used to calculate the slope of Sac atliw= 0, and, 

in fact, gives incorrectly this slope except when y = 1 


(see Eqns. (4.38)-(4.40) below). Similar remarks apply 


TOA 2s atc =" 1. 


Gy) Sonar him Cree lb 


TfO,ODcain 3s and its slope near c = 0 andce=H1, 


CC 
det us go back to Eqn. (4.13) with K # 0 and rewrite it 


in the form, using (4.12), 
Vv +V— 
(seams) ae c- vn3) = Kn,[1- (utv-1)n3]" - et 42 oN 


Consider first the case c << 1. Then noting that Nn 3 is 
necessarily less than or equal to c/u, we have, to lead- 


ing order in small quantities, 


oe 
(c - un 3) = Kn, 
or 


yo + K'y =oc = 0, (4.30) 


where 


y= (ung)7Y , xt = (yyy l/4 (4.31) 
For p = 1, (4.30) gives 
Ye Bante /dl ik); fora ae (4.32) 


For uw 2 2, we can neglect, for sufficiently small c, the 
first term in (4.30) compared to the second term and 
hence 


Vir ic) K* or ite ae cl y/K , limurce ae (4733) 
(4.33) is clearly valid if (c/k')# <A OL 
o<< (K/py te Tl? (4.34) 


Noting ‘that (4.32) and (4.33) are correct*to 


the power of c exhibited in them, and using Eqns. (4.12), 


we can readily expand the expression (4.19) for (s.8)am 
znApowers or ¢c. > "One finds, for wH*L;, 

mit coi oe 

Sac =) Cie te 2 CLR OCC) ee the (4.35) 


and*for yy 2>°2 


Siac ra i-x WG ag Sree : (4.36) 


For comparison 
Ves) sescl 


2 
= te SP 5S OnOOkS 4.37 
(Sac) Gee ct C ( ) 


2 
From (4.35) and (4.36), we have, to orderc , 
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2 
+ 0(c?) ; u 


s. 2 2v 
2 2 
Soo =O 7c + (2/K)c*+0(c7), pH2, (4.39) 
CaO ie 2 + O( 2 
cc =e Cc Ca Mas Leis ae, (4.40) 


The expressions for Sac and Sac near c = 1 are 
Simply obtained from Eqns. (4.35)-(4.40) by replacing 


everywnere Coby ti-=c, Aeby ) ‘and Vv <by Sy. 


4.4.2 Discussion and Numerical Results 


Eqns. (4.35)-(4.40) show firstly that in every 


case as c>0, S Clo as it Should. ~ Secondly, as 7c 


Ce 
increases from zero (but still c << l), Sac increases 
and lies below the ideal value if wp = 1 and above it 


Pierre. With turther, increase Ob .¢, Han. (4.27 jesiows 
that Sac has a maximum at some concentration, say, oF 
(Cc, <w/ytv), and then drops to zero at ¢c = pi/ityv. Sac 
has a discontinuous slope at this point, since as c 
increases paecnee Egn. (4.28) takes over and Sac rises 
again to a maximum, at the concentration, Coe say, and 
then decreases to zero, the limiting behaviour near 
c = 1 being governed by whether v = 1 or v 2 2 in 
accordance with equations similar to (4.35)-(4.40). 

It is interesting to remark here that the fact 
that (for K>+0), Sac = 0 at c= L/utv and has a V-shaped 


variation with concentration near c = uU/utv is similar 
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to the behaviour of the concentration fluctuations in 
BupeELlecely ordered solid alloy near the appropriate 
stoichiometric composition [Krishnan and Bhatia (1944), 
bootie cand Phornronm 3(19/0)i) aekFor the dlattergateathe 
SLOICNLOMEErIC pcOMpPOSHEION itself, of course, Sac (@) = 0 
for all wave vectors q except for the spikes which 
cCornrespond to Bragg reflections due to the formation 
Gpmencesuperlatei ce. selnid Compound forming wmmoLeen 
solution at the chemical composition, Sac (a) may then 


be expected to be much smaller than Sp Lom Sina LiaG, 


have a peak higher than aes at some higher q and then 


oscillate about oe with progressively damped amplitude 


ce 
with the sum rule for Sac (a) Given eDysbOn mec LD) a. 
Returningsito Bons ..8(4.27)- (4.28) 7 we may deduce 


from them that the positions Cc} and C5 of the two maxima 


in S., and their heights are given by 
CEP SGA) 2 ONY, Crane Oy MENA ey tea) 
Sag (Eq) = (u/y) 7£ ly) (if (y)) Cyt (y) (1-y) 1, (4.42) 
Sha(hy, = ALVA a an nae (4.43) 


where y = utv and 


2 % 
» 2y == (yo = y+ 1) 


EA) = Mma ° (4.44) 


We note that the heights of the two peaks are in general 


unequal being in the ratio (y/v) 2. Table 4.1 gives 
values of Ch, Co and the heights of the peaks for some 
values of u and v. 

For a non-zero K, but K << 1 the positions and 
heights of the peaks would still be roughly given by 
Eqns. (4.41)-(4.42), while the V-shaped behaviour near 
c= uU/utv will get rounded into a minimum. With further 
increase in K, the two maxima will gradually merge into 
each other in accordance with the limiting behaviour 


ch 


Sac =c(l-c) for K > @. These features are illustrated 


in Figs. 4.1-4.4, where we give the results of numerical 


Calculations for Sac for different values of K for the 
Caseca,)) o— (1) e372), (orl) and 27 i). 

In Fig. 4.2, we also give the experimental 
values of Soc for Mg-Bi system which is known to form 
the compound Mg Bi, in the solid state. These were 
obtained from the data on the activities tabulated in 
Hultgren et. al. (1963), by noting that if 1S is the 


activity coefficient of metal A (concentration c) in 


the binary solution, then’ 


+ Eqn. (4.45) follows immediately from (2.41) if we 
note that the activities ay and a5 in it are related 

to Yy and Yp by a, (=a,) =c¢ Ya and a, (=a,) = (1l-c) y,- 
We mention that, in certain ranges of concentration, a 
precise determination of the slopes of £n a or &n y 
versus c curves (and hence Sac! requires more data than 


is presently available. 
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Table 4.1 


The positions lon and Srp) and heights of the 


peaks in Soc for the limiting case K>+0O for some pairs 


of u,v. For comparison the value of fae = c(l-c) at 


each of the concentration is also given. 
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The concentration fluctuations Scc as a 
function of concentration un a binary 
mixture forming chemical associations 

of the AB type (u=v=1). The different 
curves are for different values of the 
reaction constant K as marked; the curve 
K=o represents Sco = Shes 6 oe) 
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Scc as a function of concentration c of 
A atoms in a binary mixture forming 
chemical: associations “Of the AB, type 


(u=3, v=2). —, theoretical for different 
values of the reaction constant K; : 
experimental from data of Hultgren et.al. 
(1963) [see also Ichikawa and Thompson 
(1973)] for the Mg-Bi system with c= Sug’ 
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Sceras#aleaunclionjoftyconcentrationge of 

A atoms ina binary mixture forming 
chemical associations of the A3B, type 
(u=3, v= 1). —, theoretical for different 
values of the reaction constant K; — —, 
experimental from data of Hultgren et,al. 
(1963) for the Ag-Al system; ---, experi- 
mental from data of Wilder and Elliot 


(1960) for the Ag-Al system with NS eat 
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Sec as a function of concentration c of 

A atoms in a binary mixture forming 
chemical associations of the A5B, type 
(u=2, v=1l). —, theoretical for different 
values of the reaction constant K; — —, 
experimental from data of Hultgren et,al.. 
(1963) for the Hg-K system with c= Cg: 
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- c(l-c) a c(l-c) 
BcCe genset Te | 3 in, (4.45) 
et, eee 
T,P Lae 


The data on free energy of mixing for this alloy at the 
temperature concerned suggests that Karrore Lepg7 4 | We 
fae that the experimental Sac has a minimum and two 
Maxima at approximately the right concentrations and 
their heights are of the right orders of magnitude. 
Fig. 4.3 similarly shows the experimental 
values of Sac for the Ag-Al system, which in the solid 
state forms the compound Ag3Al. These were obtained 
from the data given by Hultgren et. al. (1963) and also 
that given by Wilder and Elliott (1960). For this 
system the free energy data suggests an (ae and we 
Ag* z is very shallow compared 
LOrthat inawine Mg-bi systemrat Sug = 5 as expected. 


see that the minimum at c 


Finally, “an Fig. 4.4 (i=) 2,7) ve= 1), 7 the? ex- 
perimental values of Sac for the Hg-K system are given. 


Here Oe ere es We see that the height and the 


position of the first peak, and the low values of Sac 


near C,, == are fairly well reproduced by the theore- 
g 


tical curves, but the experimental curve shows no second 


Maximum. The experimental Sac for Hg-Na is Similar for 


Cuq < > but the height of the first peak is substan- 
g * 


tially less than that for the Hg-K system. These 
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discrepancies are probably at least partly due to the 
fact that these solutions, in the solid state, form 
not only the compound Hg.B (B= K or Na) but also HgB, 
Hg ,B and HgB, so that several pairs of {y,v) should 
be simultaneously taken into consideration (see 
Bons *(4.10):)°. 

Finally, as a contrast to the behaviour of 
Soc depicted in Figs. 4.1-4.4, we may recall the 
variation of Sac with concentration for the regular 


' ; Tt : 
(in the zeroth approximation) or conformal solutions. 


For these S is given by (3.34), namely, 
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te c(I-c) 
Ee Socio Gna oer 
RI 


so that it is symmetric about c = er has a maximum at 


c= > and dies, forall c,; either above sai = C=C) 


id 
or below Sacre 


energy w is positive or negative. 


depending upon whether the interchange 


t Parenthetically we mention that, as might be expected, 
the regular solution theory in the higher (quasi- 
chemical, see 83.4) approximation, gives for Soc 
essentially the same results as depicted in Fig. 4.1 - 


which is for the case ut = v = 1 and hence is symmetric 


walt 
about c = x: 


(53) 
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From the above discussion, it is apparent that 
the simplified version of the model (on which the 
calculations are presented here) has reasonable appli- 
cability to at least some of the compound forming 
solutions, although in general other factors (more than 
one type of AB and deviation from the assumption of 
ideality, see Eqn. (4.10)) have to be considered for 


obtaining quantitative agreement. 


475 Variation of Partial Structure Factors with 


Concentration 


Knowing S it is relatively a simple matter 


Cer 

to infer the concentration dependence of the partial 

structure factors ais or of the other two number- 

concentration eeructure ‘factors SxN and Suc (see §3.2). 
For this purpose we require G@vandews 26 10L these, 

as mentioned earlier, § is usually small (10.03) and 

in the absence of experimental data may, to a first 

approximation, be assumed to vary linearly with c 

between the two pure metals. 6 may be inferred from 

the data on volume as a function of concentration, which, 

if the simplified version of the model is valid, should 

be given by (4.16), in which we need only to fix the 

2°) 


fe) 
parameter v (and, of course, Vt ) and V from the 


observed data. In this connection we should note that 


according to (4.16), the excess volume AV, Pe.egthe 
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actual volume of the mixture minus the volume obtained 
by linear interpolation between the values for the 


two pure liquids, is given by 
Av=v-evy? - (1-c)vi0) = -nay . (4.47) 


Hence |AV| is maximum when nN is Maximum; that is at? the 
Concentration’ c = W/irv. For many cases” AV/V ts at most 
a few percent and we have not been able to verify this 
conclusion for lack of sufficient data. (In these 
cases, 6 calculated by assuming the linear dependence 

of volume on c forms a reaSonable approximation, cf. 
Fig. 3.8). For the amalgams discussed earlier, HgK and 


HgNa, the observed maximum in |AV| (0.24v and 0.18V 


respectively) occurs near c = = [Degenkolbe and Sauerwald 


2 
(1952), Vanstone (1911)]. We have already mentioned the 


reasons why the simplified version of the model taking 

Nee Z2and ) —"Lersunot, quate. adequate fOr Caeser systems. 
When the tendency to form chemical associations 

A BY Ys strong (Kk <<91),) the varracion of as4 wa Che con= 

centration depends markedly on the values of yu and v. 

To illustrate this feature and to see how this variation 

is affected by 6, we depict re versus concentration 


curves for a few cases in Figs. 4.5-4.7. For simpli- 


city we have throughout set 6 = 0 and drawn all the 
curves for K = TGR, Retwill “be clear fren, tie inspec= 
Emon OGst bos. 4.1454 .2OrD Sacr 2 smaller value of K 


ed. 


i | 
any I rh 
¥ : ih i : 
t i 
aay. RTT fr na a vi DY "4 
‘ ! ives i SNAG ¥ 
, wk j ' i 
' i pr i 
i Mie very | 
Wale | 
; ay ' f 
~ Wa 7 
po a 
f } {id yy! 7 
i 
i . a \ \ ‘aga. | | j 
4 ; 
“ bedi , 
i 
ah 
i ; { 
i ; Me Dae 
Cl ye 
; j 1) 
i we 
Ais \y 
Ne 5) my, } 
Wie. ca, ‘y wy hit) o 
\ | i, 7 i Ab ue ‘ is 
4 rs ies i i, f pid Mg i 
ih yd Lee MN be p cy eae the 307 iinet 
f i : i 4, e clave { val iy re yey 
y. 4 ‘ i # y Via ol 
{ ’ Le a ye ee ee ey My ee) ghee 8 5 q 
Pr ian’ ! i ,e ij we a] : is) ye ie A , . se ae wi 1% i a asst oe 
Pe, ‘ial J Nias ; 
: ‘4 


Pi Teed pil te ath Aye aay wer bar ih Aida yi. ih Ws) 

. hh RL ae Ce ks ayy! ' 

ni ee 
; d 1 ‘)) Sar ‘ ¥ 


a m 
4) dei Baie) "ane ees are A it iil Ni nat sat | 


ii Fl , 1 i 
] i eee ekg A f cn 


é ih iv j | i 
\ | r i - i j A ‘iy, i ; ) ‘ t on 
i f fe j 
wv SUN 0 , - a ae 
a eae ST Ue Ree em a 
vue ® M00, ta iy ea ae: 5 a iP oe 
| ia ; Nu AS i 
n ey 0, a Pe a iy a ee ' 
I ; ee ey ene Ney | Paid dv iih a OD ieee 8 
Ceca ees 2s aa Mei a | i U4 Vie see Oh 4 
, Rive gh ite ae | a Le yeni 
: 7 Wate Nae A w i \ | 
Wad Ph Si eee Ueki a OPE ie a 
CUM ne Ty 4 anal ” F 
iw i v ay 


bel ss BY i ay 4 eat ne ee é 3 


a Lune ‘ae 


cv in JEM 


does not affect the results qualitatively. (For K+», 


id 
Sac + Sac (=c(l-c)), and the concentration dependence 


of ass is quite different from that depicted in these 


2 
figures; for example, if 6 =-0y~then- all es = 610.03, 
andet 6 22. Or ‘ ane are as given in Figs. .3.2-3.3. 


Fig. -4,5 gives the results for the cases (u,v) = 


(1,1) and (3,1) taking 6 = 0 (equal partial molar volumes 


for the two species at all c). The curves for (u,v) = 
(1,1) are qualitatively similar to those given in Faber 
(1972). The dashed curves in Fig. 4.6 show the concen- 
tration dependence of oe for u = 3, v = 2 case taking 
6 = 0, while the smooth curves are drawn by assuming 
that the volume varies linearly between the volumes of 
the two pure constituents. The numerical values for the 
volumes were chosen appropriate to the Mg-Bi system, so 
that 6(c) varies monotonically with concentration from 
=—Q229 atecit= 0 to-—0.4 at ci= 1. Thesetvalnuesi o£ 6(c) 
should be correct to about 10% for this system since 
Av ~ 0.02 V at the composition Mg3Bi,. 

bast lv. Figuu4.7egi ves. the resulitel tor je 2, 


1 case. The curves with long dashes are again for 


Vv 


6 0, while the smooth curves are drawn using me) 


obtained from the measured data on volumes for the K-Hg 


system [Dagenkolbe and Sauerwald CL952)51 4 C6tc) first 


decreases from 6(0) =-1.3 to 6(%) = -1.5 and then 
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Theppariialescructure lactOlrsva a Basma 
function Of COnNCeENntCration Cc Of Avatoms 
in a binary mixture forming chemical 
associations of the A,B, type. — —, 
w= v=1; —,u=3, v=1. Both sets of 
curves are for 6=0 (equal partial molar 
volumes for the two species). 
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0.5 eric 
i , 


The partial structure factors a; as a 
function of concentration c of Avavoms 

in a binary mixture forming chemical 
associations of the type A3B) (u=3,v=2). 
— —, 6=0; —, 4, as explained in the 
textaj thor valuesrof appinearac= Osand of 
aggunearyc =, noticshown ingthe diagnam, 
see text following eqn. (4.56). 
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The partial structure factors aj. aS a 
function of concentration < of Asatoms 
in a binary mixture forming chemical 
associations of the type AjB) (u=2,v=1). 
, theoretical with 6=0; —, theo- 
retical with 6 as explained in the text; 
---, experimental from McAlister and 
Turner (1972) for the Hg-K system with 
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increases monotonically to 6(1) = -0.7). For comparison 
the values of 234 obtained by McAlister and Turner (1972) 
from the measured thermodynamic data are also given. 
Considering the limitations of the model discussed 
earlier and the uncertainty in determining Sac from 

the measured activity data, the agreement between the 
two is not unsatisfactory. 

The results presented above need two further 
comments. First, we observe that in Figs. 4.6 and 4.7 
the dashed (6 = 0) and smooth (6 # 0) curves for each 
ai5 come very close to each other near the appropriate 
chemical composition. This is not surprising if we 
Sreter to Eqn. (2.36) for as4 and consider the limiting 
case K>0. Then since now Sac = 0 at the stoichiometric 
composition c(=c,)) =u/utv, the partial structure factors 


become independent of 6 at this composition, and are 


given by (with 6 = 0) 


CaP es = Ee 
=-F and aj,=tl, (c “ery 3 (4.48) 


The values of ai ee OS = given in Figs. 4.6 and 
3 Z 
4.7 are close to these values at c = ¢ and c = z 


respectively. 


Second, it is of interest to give Sexplacit 
expressions for the values of aes atc. ="0 sand c=. 
If we substitute the expansions (4.38)-(4.40) for Sac 


into the Eqn. (2.36) and take the limit c > OC lorece a 
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as the case may be), we obtain ta cs Des etc.) 
1) 


ay tetes . a0 = 9° for all u and vy, (4.49) 

255 =i6P eGo nae = gt 4 5h £forvall eandaw, (4.50) 

ay) = 0° - 26° -2v(1+K) >, for p = 1 : (4.51) 
=e9° — 26° + 2/K , for p= 2 : (4.52) 
=@° - 26° Aeon ise 3 s (4.53) 

and 

ayy = 0) + 26'-2u(1+K) 1, for v= 1 : (4.54) 
= 69> + 26) + 2/k 4 ay SD ; (4.55) 
Sy ee ae 5 aes b “ase 


Thus, for example, in Fig. 4.6 (yu = 3, v = 2), the dashed 
curve for ayy starts from zero and steeply rises with 
small increase in the concentration, has a large maximum 
(N33) "not “shown in’ the diagram) “at c= 0.01, and then 
follows the curve as in the diagram. In contrast anor 
near c = 1, continues to rise to the value aoe ea 2000 
in accordance with (4.55). As shown by Eqns. (3.10) and 


(ool and Ans depend on the curvature of S,,, i.e. 


ita 
sgSeesdce, at c = 0 and c= 1 respectively, so that 


their determination from experimental Sac TS ca Cult. 


The limits given in Eqns. (4.51) to (4.56) may therefore 


serve as useful guidelines where the model is applicable. 
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CHAPTER 5 


CONCENTRATION FLUCTUATIONS AND THERMODYNAMIC PROPERTIES 


OF SOME COMPOUND FORMING BINARY MOLTEN SYSTEMS 


5.1 Introduction 


In the previous chapter we discussed some 
general characteristics of the concentration fluctuations 
Sac and other structure factors in the compound forming 
binary molten systems on the basis of the so-called 
chemical approach. This assumes that if the mixture of 
A and B atoms forms, in the solid state, a compound at 
a single chemical composition’ specified by AB, Cate Shea 
small integers), then in the liquid state there exist, 
at a given temperature and pressure, certain numbers of 


A and B atoms and chemical complexes A B, in chemical 


u 
equilibrium with one another. The equilibrium numbers 
of these and hence the macroscopic thermodynamic proper- 


ties of the binary mixture depend on (a) the (free) 


t The more general case where compounds are formed at 
more than one chemical composition was also formulated 
in the chapter, but will not be pursued here. If one 
of the compounds is considerably more stable (for 
example, has much higher melting point) than the others, 
assume that one type 


one may as a first approximation 


of chemical complexes is formed. 
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energy of formation of the chemical complexes or the 
reaction constant, and on (b) the form of the chemical 
potentials assumed for A,B and A,B, that is, whether 
the ternary mixture can be regarded to behave like an 
ideal, athermal, conformal solution etc. 

The behaviour of Sac etc. for different values 
of uw and v was discussed in Chapter 4 by making the 
simplest, namely, the ideal solution assumption referred 
to in (b) above. However, this assumption is too crude 
to give, in general, quantitative agreement with experi- 
ment. In this chapter we discuss two higher approxima- 
tions and calculate on their bases the concentration 
.dependence of the free energy of mixing, the activity 
and Soc for four binary systems: Ti-Te, Mg-—B1, (Cu-Sn 
and Ag-Al. Each of these higher approximations contains 
four interaction parameters, and although these have 
physical interpretation they have to be at present de- 
termined from the thermodynamic data themselves. 

The two approximations are described in 85.2 
and the various relevant formulae are derived on the 
bases thereof in 85.3. In §5.4 the behaviour of Sac 
is discussed for two limiting cases for which the 
equilibrium equation can be solved analytically. The 


results of numerical calculations and comparison with 


experiment are presented in 85.5. This is followed by 


a brief discussion of the results and some concluding 


remarks. 
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5.2 Basic Formulae and Approximations 
eee eee PPL OAIMActOnsS 


Let the binary solution contain in all Ny=wWe 
and Np = W(1l-c) gm moles of A and B atoms respectively, 
c being the atomic fraction of A atoms. We assume that 
only one type of chemical complexes Bes (u,v small 
integers) are formed. Then if there are in the solution 


nj, gm moles of A atoms, no of B atoms and n, gm moles..of 


Bree we have from conservation of atoms 


nj=we- uns , no = W(1-c) - vn 
and 


n=nj,+n,+n,=WV- (utv-1)n, : (Sel) 


The free energy of mixing for the binary A~-B alloy can 


obviously be written as 


al Soy Led Py (o) 
Gy = G -~WcG, Wl c)G. 
suier eewins Gs (5372) 
with 
g = nee) + yee - ne (Si53)) 
ate aes ee ets fae (5.4) 


where glo) i = 1,2,3, is the chemical potemtial for 

ah 
the pure species i in the solution. The equilibrium 
value of n3 at a given concentration, pressure and tem- 


perature is given by the condition that G and hence 
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(9G e) = Q-. 
Reg ts eats Te (5.5) 


We remark that in Chapter 4 we had formulated the problem 
somewhat differently (using the language of chemistry 
rather than physics!). The equivalence of the two 
formulations is demonstrated and their relative conven- 
ience commented upon in Appendix 1. 

Erotss2)t ehevrfirst ccexrm (-n3g) represents the 
lowering of the (free) energy due to the formation of 
the chemical complexes. The second term is the free 
energy of mixing of a ternary mixture of fixed Ny Ny MN, 
whose constituents A,B and ee will be assumed to inter- 
act relatively weakly with one another - the strong 
bonding interaction between A and B atoms having been 
already taken care of via the formation of chemical com- 
plexes. We may thus borrow expressions of varying com- 
plexity for G' from the various theories of weakly 
interacting mixtures. 

The simplest expression for G' comes from 
assuming that the ternary mixture forms an ideal 


solution whence 


3 
G' = ‘RT bs n, gn (n,/n) ; (5.6) 


which is just (-T) times the entropy for random mixing. 


The discussion of concentration fluctuations on the 
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basis of (5.6) and (5.2) was given in Chapter 4. Eqn. 
(5.6) is valid when (a) the effects of differences in 
sizes of the various constituents in the mixture can 
be ignored and (b) the differences (defined below) 
in the interactions between the different species are 
zero. If we retain the assumption (a) and treat the 
interactions whi to be small, then we have conditions 
under which the theory of regular solutions in the 
zeroth approximation or the conformal solution appro- 
ximation is valid. G' then is just a straightforward 
generalisation of (3.24) to a ternary mixture and is 


given by [see, for example, Longuet-Higgins (1951) ] 


fe y 
GuaeRTS) n, gn(n,;/n)+ ) } (nyns/n)o,. (527) 
aca 
where 3,3) = 12,3, and Wi4 (=0, for 2) =) j) sare toe anter— 


action energies defined in tie usual way, for example, 
2015 = 2Wag 7 Yana Op etc. For brevity, we shall refer 
to (5.7) as the conformal solution approximation or 
simply approximation (a). 

The effects of differences in sizes between 
A} BoandsA B are more difficult to take into account. 
A simple approximate expression for G', in similar vein 


fy (527)-418: that due to Flory (1942, 1953), well known 


in polymer physics, namely, 


GieeRT Jin, (inyey hd 2 ns 5X G5 : (5.8) 
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where 5 is the concentration by volume of the ith 
species in the mixture and Xi (20015 15> Sy eee 
interaction energies (similar to oe in, Eqn. (5.7)) 
between the different species, and are (RT) times the 
X45 eefinedjingFtory.,(1953)>. The) fi¢cst+termvin: (5.8) 
is the expression, in Flory's approximation (see also 
below), for (-T) times the entropy of mixing three 
species whose molecules differ appreciably in volume 
from one another. If these volumes are all equal, then 
5 = n,/n and (5.8) becomes identical with (5.7). We 
make here a further simplification of (5.8) by assuming 
that the volume per atom of A and B atoms is nearly 
the same, say v, and the volume of AUB se (ery) vi. 
[This assumption is essentially similar to that often 
made in quasi-lattice models of monomer-polymer solutions 


and is of course only approximately true!; see below]. 


Then 

Pe /N ois O50 = NoLVe oy) | a ae OEE Sod) 
With (5.9), expression (5.8) becomes 

G' = RT[n, fn(n,/V) + nz kn (no/WV) + ng Pn ((utv)n3/N’)] 


(DLO) 
peer ey ei i ie 
ks} 
Jaren nnn EEE 
uy For an account of the quasi-lattice models, see, for 


example, Guggenheim (1952) or Prigogine (1957). 
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where we have set yee oe a (utv) Xj 3¢ Yo = (u+V) X53- 
We shall refer to (Solos for brevity, as Flory's 
approximation or approximation (b). We note, as is 
well known, that the combinatorial term [the term 
proportional to RT] in (5.10) follows also from the 
quasi-lattice models of monomer-polymer solutions if 
the coordination number z is taken to be infinite’, 

In the following we shall use the approximations 
(a) and (b) above, in conjunction with (5.2), to calcu- 
late the thermodynamic properties of some compound form- 
ing systems. Each of these approximations contains 


four parameters, g and w.. 


+5 (or Vis). There exist, of 


3 

course, hagher approximations to G* than (5.27); {5238) 

or (5.10) in the literature, but these are involved 

and contain, in general, additional parameters 
[coordination number, surface to volume ratio for polymer 


molecules (chemical complexes) etc.] and will not be 


pursued here. 


t See, for example, Guggenheim (1952), Eqns. (10.10.4) 
ana (10.10.11). ° Note: that in these equations 7far our 
problem, ae te r= nb r3= utv. Note also that N; of 
Guggenheim are n; in our notation, and the combinatorial 
termain G’ 2s just, —T } n,AyS- In the limit z>~%, 

G37 uty and the (10.10.4) gives the combinatorial term 


in (5.10). In practice, the difference between the 


latter and (10.10.4) is quite small for z> 8 or 10. 
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It may at first sight seem surprising that 
we should consider the approximation (a) at all, since 
even if the sizes of A and B atoms are equal, the size 
of A, B,, will not be that of A or B. However, we should 
brecall firstly from §3.4, that the conformal {or 
regular) solution theories aie generally considered 
to be reasonable for volumes of the different species 
differing by up to 100%. Secondly, the bond length 
between A and A, B and B or A and B may be different 
(less!) when both atoms of the pair are in the chemical 
complex than “spi: B= or both of the atoms do not belong 
to the chemical complex. Hence, since approximation (b) 
is also only an approximation to a truly complex situa- 
tion, and uw and v are small integers in ook problem, 
it was thought desirable to make calculations on both 
the approximations (a) and (b) above. We comment later 
on the extent to which the two approximations agree 


with experiment. 


5.3 Expressions for Various Thermodynamic Quantities 
Exp ess 1 Os Se 


530i) sApproximation (a) 


Combining Eqn. (5.7) with (5.2), the free energy 


of mixing Cu for the compound forming binary alloy is 


Gy=- n3gt RT ee n, &n(n,/n) + bl (njnj/n)ws, + (5-11) 
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The equilibrium number of chemical complexes is given 
by the condition (5.5). This gives, on using, (5.1) 


and after some rearrangement of terms, the equation 


)Ke (5212) 


where 


K = exp[-g/RT] (SiS) 


and 


Y= (w) 5/RT) { (utv-1) (n)n,/n*) - u(n,/n) - v(n,/n) } 
+ (wy 5/RT) { (wtv-1) (nyn3/n*) - u(nj/n) + (n/n) } 


* (Ww 4/RT) { (wtv-1) (njn,/n*) - v(n3/n) + (n/n) }. (Seis) 


On eliminating Ny, Ny and enyin).(5.l2)andy 45) B44), from 
(5.1), Eqn. (5.12) is seen to be an equation in the 


Single unknown n When the solution of (5.12) for 


3° 
Nor and the corresponding values of Nj, Ny and n are 
substituted into (5.11) one obtains the equilibrium 
free energy of mixing Gur for comparison with experiment. 
In the following, to avoid cumbersome notation, we 
shall continue to denote the equilibrium values of N3, 
Nj, No, etc. by the same symbols. We note that the 


differentiation of Gur with respect to a variable say T, 


at equilibrium can be written as 
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(9G,./9T) = (3G /3n.) (dn,/9T) i 
M P,cw MW 3° 0, PCN cl P,c,W 
+ (0G,/dT) 
M P,C,n3,N 
= (3G /9T) aye is 
u/ P,C nN 


by virtue of the equilibrium condition (5.5). Similar 
testlts apply for differentiations with respect to c or 
sed 

Using (5.11) and (5.15), the expression for 


the heat of mixing H,, is 


M 
Sh a eRe Se) (S2i6a) 
ag wee OW is 
==_— =- aot ed Rar td 5 SLED 
n,(g r(3) | +) L (")) 4,5 -Tp Le ( ) 
i<j 
and the entropy of mixing Su = (Hy 7 Gy) ZT: 


The concentration fluctuations Sac are given 


by (2.39), namely, 


2 2 
wea BY; 
Soc = WRT/ (9 Gy/ oe ) .P,N? ( ) 


Hence differentiating (5.11) twice with respect to c 
and using (5.1) and (5.12), one may obtain after some 
algebra 


Ss 
S Seen ee (5.18) 


ce 
iL + & Sac 


Here 
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3 ~1 


Sag =W a (nt)?/n, - (n')2/n : (5.19) 
Leese Ce ns ’ 
& = aa: r) wi A (5.20) 


where a prime denotes differentiation with respect to c. 
The expression for ny is obtained by straightforward 
ditferentiation-of (5.12) using (5-1) and moting that 


nj = WV - n3 ete. Ins (52.20) 
(n./n) > = min = nen} /ue é 
a i 


Finally, as mentioned previously, experimen- 


tally Ss is determined from the measured activities. 


EC 
a fe an denotes the activity of the component A in the 


mixture, one has from (2.41) remembering that ay= ays 


ws as ' 
Sac ae (1 c)a,/ay e (S221) 
In terms of Gur an tsigiven by Eqn. (2740)i, swe". 
dG dG 
M 1 M 
RT &2n a,= - i{< + (l-c) (x— te Ot) 
A oN, 1,P,N N| M dc T,P,N 


B 
With Cu given by (5.11), fan. (522)) gives, onusing 


(S275)p and (>)..9), 

n nN n nN. nN. 

1 i 3 aa 1 a : 
en ays in() + ppl Mat G@Mia7 bd BD qba,5| 


(S223) 


Using (5.23), the expressions (5.21) and (5.18) for 


Scc may be verified to be equivalent to one another. 
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The expression (5.18) has certain advantages in deducing 
the limiting behaviours of Sac: considered in §5.4, from 


theory. 


5.3.2 Approximation (b) 


In this approximation G' is given by (5.10). 


Hence using (5.2), we obtain! 
Gy = = n3g +RT(n, gn WwW No Qn Fie 3 gn ad 
+ ) ) (nen MW Jv, : (5.24) 
7] 


Using (5.5) in (5.24), the equilibrium value of n is 
given by the equation 

x An expression essentially similar to (5.24) for Gy 
has been recently derived in another manner by 
Takeuchi et. al. (1973) for the case where one of u,v 
is.unity. Takeuchi-et. al. .use this, toscalculate 

Nn, as a function of concentration for Cu-Sn and Cu-In. 
The interaction parameters in this work were deter- 
mined from the observed values of Au neglecting in 
the expression for Au the possible temperature depen- 
dence of Vise No comparison of other thermodynamic 
quantities is made in this work. (A calculation made 


by us for Cu-Sn with their parameters gives substan- 


tial disagreement with experiment for Sac: See also 


7552 
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nin, = (nw ekene (529259) 
where 
K' = K(utv)exp[- (yutv-1)] (3'2'6)) 
Z= (NRT) 


[ (nj -un3) vy 3+ (ng-vn3)v53-(ungtvn,)v) 9] 
(5er25/,) 


and K is given by (5.13). Next remembering that the 


heac-of maxing and activity are related to Gu bye Abii 6.a) 
andre 5-i22)% one’ obtains: fore this: case 
ag ; ; n,n. ee 
H, = -n [g- T (—=) | + Board v5 Ta (55.28) 
M 3 di) p i<j N 2 OD 
n 
com A Ne 1? 1h 
Qn a, = £n aay aa WN WRE (N3V)3 + N5Vj 9) 
- = ye) n.Navis . (529) 
WORT i<j J +4 
Pinay ~susting, (os45/) tia et SecA) mono melds 
that Sac may. again formally be written in the form 
Coeeuc 7 bike wlth Sac and a in dt neweoiven py 
toy?) 
Ss aie (Gaus Ane 
ee q=4) x 1 
Zz ro 
= 2 ACR as ae) njnivig + (53 ON) 
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5.4 Expressions for Sac for Two Limiting Cases 
(22 IES Da Rr an ccs Seth ME ER TD ED 


Of the thermodynamic functions Gur an and Sac: 


the concentration dependence of Sac exhibits the most 
interesting and varied features (see Figs. 5.1-5.3). 
To have a qualitative understanding of these features 
it is instructive to consider two limiting cases for 


which analytic expressions for S may be obtained. 


ce 
The first limiting case refers to situations 
where there is a very strong tendency to form chemical 
complexes, that as, og/RT >> 1 or \K << 1s... We “therefore 
consider the limit K*0. Considering first the appro- 
ximation (a), we see that in the equilibrium Eqn. (5.12), 
K> 0 implies also K exp[Y] 70, since (oS 57 are by 
assumption small finite numbers. Then, if nN 7 0, 


Eqn. (5.11) implies that either nj? 0 or no* 0. One 


then has, using, (5.2), 
n> Wc/u for 0< com Loi Vy (oar) 
ee NV (1-c) /v £or U7 (OY) es = Se 4 (5232) 


ine tne Vamrt (5.31) n,7 0, and one may eiow also that 


t To see this, let nz = (Wec/u) - An. Substitute this 
ini(cb2). <Then using -(S.1) , rememberimgn that hss ak, 
and retaining lowest order terms in An3, one finds that 
An, and hence n, =pAdn, are proportional to a. positive po- 
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wer of K. Hence as K~0O, Nyy and (n,) /n, tend to. Zero, 
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n,7 0 and (n})?/n,> Oe eUsing this! resultwands(5e3)) 
andeigns.| (521) inthe expressions (5419) sand .(5.20) 
for Sy, uke &, one obtains straightforwardly for 
Ome ee U/iCuty)) 

oar eee / W)C Cu eeu (eva 4 ).C] 


B'9) = -2 (0 4/u7RT) (Wn) 3, n=NE1-c(ytv-1) /ul. (5.33) 


The index 'a' signifies that these expressions are for 


the approximation (a). Be then is of course 
(a) (a) EVE 
Sac = Sac [l + & Sac ] : (52/34) 


Simivarly trom (5.32)),.f£or the concentration. range 


W/E) << Co <1, 


ee cli y= (lee lve (lecito en 


G?) = -2(w,4/v RTI (W/n)? , n=WU1= (Ae) (utv-1) /v1. 
(558359) 
The solutions: (5. 31) sand! (532)Pand hencesthc 
corresponding expressions for Sac and) 2 are: not avalid 
in the limits c> 0 and c+l1 respectively. For these 
cases n, itself tends to zero and the equilibrium 
bone. (5.12) srequiresamore careful analysis. By follow- 
ing steps similar to those given in 84.4 for the case 
Way = 0, one may show that for c << 1, the equilibrium 
Fane. (5512) has. the solution 
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2 
nz = We COMO (C1) Pe PLO eal (5. 36a) 
n 3 = We"/K, a3 One) ASTLor oe (5563 '61b) 
where 
ee a ta -] 
K, = K expl (w,., WW) 5)/RT] ier ie tiie (1+K.) C5 3m) 


livemexponent) 1m. 5.37) yisy just. theivalue or Y Eeqna(>.14)| 


arco = 0. The solution (5.366) is valid iwhen both (see 


84.4), 


Cie | and Cu < CON Ape oe (ye She!) 


As may be verified, the solutions (5.36) are sufficient 


to determine oie! GoRorder co? and one obtains after 
some algebra 
gan! ae 5 - e 2va + a (W134 (1-a) -w5 30-wW, 5 (1-a) ) J 

+ O(c?) for u=l, (5.39) 
and 
g (a) = ee N SR Ob ar )-2w /RT] + O(c7), £oGg, Us 42 
€e a 12 


(5.40a) 


= Bee Gi - 20, o/RT] + O(c), EOI eo ae. (5240) 


The expansions (5.40) are; valid when both the inequali- 
ties in (5.38) are satisfied. *,To obtain similar -expan- 
sions near c = l, replace in Eqns. (5.36)-(5.40), c by 
l-c and make the interchange u Z v and W532 W413: 


Eqns. (5.39)-(5.40) demonstrate that for sutfivcientiv 
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small c, Sac = c in all cases. The expansions are 
carried to order ce here, since as discussed earlier, 
the value at c = 0 of one of the Faber-Ziman partial 
structure factors depends on fae eer / dc alee 

Similar remarks with appropriate changes apply to the 
region: C.1.; 


We note that in the special case where all 


Wis = 0, so that G' is given by (5.6), the various 
Expressions (50a 3)45 35008 (5009), and (5.4ny for eae 


reduce to those given in Chapter 4, as they should. 


In Flory's approximation (approximation (b)), 


LOL Kk: 2:0), n3 are again’ given. by? (5.31))and 45. 42).. 
For Sacr By etc. cone then obtains, using, (a.50) > 

(Di pe Cd us (uty) Cc) (1) Seg eit) 5 4 
Scc ~ ut (uty) (utv-lL)c Tacos tie u2RT rs iay (5.41) 


Lore Ounce) <au/ (pt) ,, cana 


Oe Wile) oy Waa) (b)__ 2(utv) 
Sis = Sco geran Werraren hy Glcs) Gane <) SoHE lc be (5.42) 


for p/(utv) < c < 1. Finally the low c expansions for 


this case are 


s(P) = exo? { (1+v8) 7+ (2/RT) [vy 9 (1-8) (14VB) + (1-8) BV, 


SPUivvelevaal) + Ole), for ved 9) (oeaa) 


e-c? [1-(3/K,) -2v, 9/RT] + O(c") >. Hee dns teitones Ean 


c-c7 [1-2v, 5/RT! Paice as fay SA) 
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where 


K, = K'exp[(vy3-uv,5)/RT] , 8 = (1+K,) 7? : (5.46) 


The expressions (5.44)-(5.45) are valid only when both 
G<< lvanditc << Ay gait & The expansions near 
c¢ = 1 are obtained from (5.43)-(5.46) by interchanging 


c i (l-c), wz Vv. and. v 


ed ali 

The second limiting case refers to situations 
where there is only a very weak or no tendency to form 
enuemecal ~complexes, “that is, (g/RT) << OtopeK «>> i erAs 
Kasea, Sain sanspection of the equilibrium: Eqni:. (5.12)mor 
(5.25) shows that in both the approximations (a) and 


(by naka Kr, Hence as K > », nj, n4 and ft) Sm. 


3 
all tend to zero. Then since n)> Wc, no > WV (1l-c) and 
n->wW, we have for approximation (a), from Eqns. (5.18)- 


poe 20)') 


c(l-c) 
= K - 5.47 
Scc 1-2 W2/RT) c (1-c) : ( ) 
The approximation (b) gives the same result with W192 
replaced by Vi2° The expression (5.47) for Sac aK 
the same as for a binary regular solution in the zeroth 
approximation, as is to be expected. 
From the various expressions given above we 
infer the following regarding the concentration depen- 
dence of S..: (1) In every case Sac starts from zero at 


CCy 
c = 0 with the slope unity and ends at zero with the 
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slope -l at c= 1. (Z)38f Ko << 10 then sinethe: kK 0 
limit Sac is again zero at the chemical composition 

cy = u/ (uty) = see \Bqns..(5.55)— (5.35) pane a(504) = 5 242)9. 
Thus Se has two maxima, one in the concentration range 
O0< cee and the other in the range ee ees Eo elne 
height and position of these maxima depend, apart from 
mene values ‘of K(<<1), 1 and v, on the magnitudes and 
Signs of Wo 3 and 013 (or V53 and V3) respectively. 


(3) In the opposite limit K>»~, § has just one 


cc 
maximum which occurs at c = sr its height being deter- 
mined by the sign and magnitude of W192 (on Vi2)> We 
note that an infinite Scar which can happen for posi- 
tive Wis (or Maa? implies phase separation and our 
treatment applies in such cases only at temperatures 
above the critical temperature for phase separation or 
far away from the critical composition. 

The results of the detailed calculations show 
that the K>+0O limit expressions for Nn and Sac form a 
good first approximation if g 2 3(yt+tv)RT. Remembering 
(5.31) and (5.32) this implies that, with Gu negative, 
1G, (e) JN > 3RT. Two of the systems (Mg-Bi and T1-Te) 
examined below belong to this class at the temperatures 
under consideration. For the other two systems 
IG (C.) /WRT | is not so large and the behaviour of 
S is intermediate between that expected from the K+ 0 


Gl Gs 
and K > ~ limits. As the temperature increases 
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K (=exp[-g/RT]) may be expected to increase, so that, 
for a given system, with increasing temperatures, Sac 
will tend towards the behaviour given by (5.47) and 
eventually towards Sac = c(l-c) appropriate? to-an 


ideal Bed utionus 


5.5 Numerical Results 


In order to apply the various equations of 
§5.3 to a specific system, we have to solve numerically 
the equilibrium Eqn. (5.12) and at the same time deter- 
Mine the parameters g (or K) and ae from the observed 


data on the free energy of mixing G SATs De ele 


M’ 


(We explicitly confine our remarks to approximation (a) - 
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the procedure followed for approximation (b) was similar.) 


: We note that the fact that the two approximations (a) 


and (b) give the same limiting behaviour for Sac at 
high temperatures is a consequence of the assumption 
that in both the approximations the atomic volumes of 
pure A and B liquids are taken to be not too different 
from one another. If this is not the case one should 
use for G' expression (5.8) and then the limiting beha- 


viour of S would be given by the expression which one 


Gi 
obtains by considering the binary mixture in Flory's 


approximation, i.e. Eqn. (3.44) of Chapter 3. 
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For two of the molten systems considered, 
namely Mg-Bi and Tl-Te which in the solid state form 
highly stable compounds Mg Bi, (u = 3, v= 2) and 
T1,Te (u = 2, v= 1), IG, (c.) WV RT | > 3! atthe ‘tem-— 
peratures of observation. For these cases, therefore, 
we took as a starting approximation that n is given 
bye (5-51) and (5.32) which! is’ appropriate for kK << 1, 
Then at the chemical composition COSTE), Gu * “N39, 
which gives a starting value for g or K. With n3 stil 
given by (5.31) and (5.32) and g as determined above, 
the expression for G,, for c <Cyr contains the single 
unknown parameter Wo3 (since n> 0). A value of Wo 3 
was thus determined from the observed data on Gur at an 
intermediate concentration between 0.1 and cu: 
Similarly W13 was determined from the observed value 
of Gu at a concentration intermediate between Cy and 
0.9. It may be seen, using expressions (5.31), (5.32) 


ard (5.50) lor mn that for systems with K << l, W412 


3" 
affects Cu significantly only for Ue t2sandect<< was 
determined by the second inequality in (5.38) [and 
simi tarly tor > 2) an -thesegion l-ee <i eoiunce 
this implies c to be less than 0.001 for Mg-Bi and 
Tl-Te systems, w could not be determined from the 


ie 


experimental data on Gus Finding that any reasonable 


choice of WI (vRT) affects the values of other para- 


meters only slightly we set in our calculations Wi = 0. 
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With the above starting choice of g, and Wi5r the 


equilibrium Eqn. (5.12) was solved for n. and the 


iS 
parameters were adjusted in the light of these values 
of n. COVODtain ia good’ £1.t° for Gus The final choice 
of the parameters is given in Table 5.1(a). 

For the other two systems considered, namely 
Ag-Al and Cu-Sn, which were assumed to form chemical 


complexes Ag Al Ci ="3, 0 = 1)" and Cu Sn Goes v= 2) 


4 
respectively, |G, /VRT| = 1.2 and 0.9 at the respective 
chemical composition. This implies smaller values of 
g/RT (larger K or weaker tendency to form chemical 
complexes) than for Mg-Bi and Tl-Te. As a starting 
approximation we took g to be -(ut+v)Gy(c.) AN. Unlike 
the case of Mg-Bi and Tl-Te, the calculations had to 
be made for several different choices of parameters 
before a good fit for Gur could be obtained. We did not 
attempt any mean square deviation tests to decide on 
the best fit so that the parameters given in Table 5.1 
are to be regarded as reasonable choices rather than 
necessarily the best possible ones. 

The parameters similarly found for approxima- 
tion (b) for the four systems are given in Table oie wl ed le) 

The results of calculation of Gy together with 


the experimental data for the four systems are plotted 


in Figs.5.1 (a-d). The variation of nN with concentra- 


tion is also shown here for each case. Figs.5.2(a-d) 
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Table- 5.1 


Interaction parameters for the systems Tl-Te, 
Mg-Bi, Ag-Al and Cu-Sn in the approximations (a) and 
(Db). 


(a): Approximation (a) 


21. Oue 


Pe252 1.0% 
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See discussion in the text regarding the choice of Wi 


and v for Tl-Te and Mg-Bi. 
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ss Ro pa Ms Ba y= Gy/W RT and n./jV as functions of con- 
centration for Tl-Te. Curve A: conformal 
solution approximation; curve B: Flory's 
approximation - see text. ® Experimental 
data for Gy taken from Nakamura and 


Shimoji (1971). : 
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Fig. 


ayia B iA 


Gy/YRT and n./fVY as functions of concen- 
tration for Mg-Bi. Theoretical curves A 
and B as explained in Fig. 5.la. e@ Experi- 
mental data for Gur taken from Hultgren et, 
al. 11963). 
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Gy/WRT and n,/V as functions of con- 
centration for Ag-Al. Theoretical 
curves A and B as explained in Fig. 5.la. 
@ Experimental data for G,, taken from 
Wilder and Elliot (1960). 
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Gy/WRT and n,/V as functions of con- 
centration for Cu-Sn. Theoretical curves 
A and B as explained in Fig. 5.la. e 
Experimental data for Gur taken from Hager 
eteals (1970). 
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Pig. S.2a. The logarithm of! the acta yinyey 7 any: 
versus concentration for Tl-Te. 
Curve A: conformal solution approxima- 
tion; curve B: Flory's approximation - 
see text. e Experimental data from 
Nakumura and Shimoji (1971). 
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Fig. 5.2b. The llogarithm of ‘the activity, In aug 
versus concentration for Mg-Bi.- i 
Theoretical curves as explained in 
Fig. 5.2a. e@ Experimental data from 
Hultgren et.al. (1963). 
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The..logarithm of the activity,.in ange 
versus concentration for Ag-Al. g 
Theoretical curves as explained in 

Fig. 5.2a. e@ Experimental data Pegs 
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In Oc, 


Fig. 5.20." The logarithm Of the activity, an Aay! 
versus concentration for Cu-Sn. 8 
Theoretical curves as explained in 
Fig. 5.2a. e@ Experimental data from 
Hager et.al. (1970). 
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The concentration fluctuation Scc versus 
concentration for Tl-Te. Solid curve A: 
conformal solution approximation; curve 

B: Flory's approximation - see text. 

— — — Experimental using data from 
Nakamura and Shimoji (1971). ----- 

Soc l=e(1-c) ] for an ideal binary solution. 
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Fig. 5.3b. The concentration fluctuation Sac versus 
concentration for Mg-Bi. Theoretical 
(solid) curves A and B as explained in 
Fig. 5.3a. — — — Experimental using 
data from Hultgren et,al. (1963). 
----- S.-~l=c(1-c)] for an ideal binary 
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Fig. 5.3c. The concentration: fluctuation’ Soc versus 
concentration for Ag-Al. Theoretcial 
(solid) curves A and B as explained in 
Fig. 5.3a. — — — Experimental using 
data from Wilder and Elliot (1960). 
----- Sacl=c(1-c)] for an ideal binary 
solution. 
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2 3d. The concentration fluctuation Sce versus 


concentration for Cu-Sn, Theoretical 
(solid) curves A and B as explained in 
Fig. 5.3a. — — — Experimental using 
data eae Hager et.al. (1970). 

----- [=c(l-c)] for an ideal See 
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compare the activities and Figs.5.3(a-d) the concen- 
tration fluctuations Sac: For Tl-Te the height of the 
second peak in the theoretical curves is very sensitive 
to small increases (1%) in 4 3/RT (or V43/8T) beyond 
the values listed in the Tables, indicating the tendency 
for phase separation. For the values of 013 (and V13)" 
given in the Table, the height of the second peak (not 
aati in the diagram) is about unity. 

It is well known in the theory of binary 
regular or conformal solutions [see, for example, 
Guggenheim (1952), Longuet-Higgins (1951)] that the 


concentration dependences of G and of the heat (Hy) 


mM! 
and entropy (Si) of mixing can in general be simul- 
taneously fitted with the experimental data only if 
the interaction parameter is assumed to be temperature 
dependent. In the same spirit in our work we have to 
assume that g and Wis (or Vas) are temperature depen- 
dent as we have done in deriving expressions (5.16) 
and =(5.28) for Hy: If the temperature dependence of 


gand «.. (or ae is ignored then Hu and Sy are 


ij 
readily calculated using the values of n3 [and con- 
sequently of ny and nN from Eqn. (5.1)] given in 
Figs. 5.1(a-d) and the values of g etc. from Table 
Sa. “Onewtinds. thatucor o1-Te, the wvalues or Hy S° 


calculated agree with experiment to within 5% over 


the whole concentration range, while the deviations 
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for Su are considerably larger. In particular at the 
chemical composition c = 2/3, the calculated Su =~ Q, 


while the observed Su is negative. In contrast, for 


Cu-Sn system both Hv and Su so calculated are consi- 


derably in error with experiment. The relative 
importance of the terms involving the temperature 


derivatives of g and Ws (or Via) in the expressions 


for Au and Su would of course vary from system to 


system, and the above comments are made here only to 
illustrate that the parameters g etc. should be deter- 
mined from the observed values of Gu rather than of Au: 


As an illustration of the extent to which 


agreement can be obtained for H,, and Sur we give plots 


M 
of Hy and Su for Tl-Te (Fig. 5.4), in which the theore- 
tical curves were calculated by taking dg/dT = -0.44R, 
and 
dw, 3/daT = =1.2R., for approximation (a) 
and 


dv, 3/dT SA eh ag dv,3/dT= -R , for approximation (bh) = 


(5.48) 


The temperature dependence of the parameters not listed 


was taken to be zero. The fact that Hu 1s concave for 


c<c and convex for c> Co is connected with the fact 
c 


that Wo 3 (v.53) and 013 (v} 3) are of opposite signs for 
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Heat and entropy of mixing for Tl-Te 
system. Theoretical curves A: conformal 
solution approximation and B: Flory's 
approximation - see text. e Experimental 
data from Nakamura and Shimoji (1971). 
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this system. Note also the signs of their temperature 


derivatives given in Eqn. (5.48). 


9-6 Discussion and Concluding Remarks 

We observe from Figs. 5.1(a-d), that both the 
approximations (a) and (b) reproduce the observed con- 
centration dependence of the free energy of mixing Gu 
remarkably well for the systems examined. The activi- 
ties and Soc which depend respectively on the first 


and second derivatives of G,, with respect to concen- 


M 
tration provide more sensitive tests for the model. 
As seen from Figs. 5.2 and 5.3, the agreement with 
experiment for these quantities also is not unsatis- 
factory. For the Tl-Te and Mg-Bi systems for which 
Wi > and Lp hake set equal to zero, the theoretical 
values of the activities a and aug can be somewhat 


TL 

adjusted in the concentration range ONsiC 3C8 by taking 
a non-zero 05 and V2 [see Eqns. (5.23) ana (5.29) e 
but this does not seem to improve the overall fit with 
experiment. 

As regards the values of the interaction 
parameters required to fit the experimental data it 
is gratifying to see that for each of the systems 


examined, the two approximations yield for the (free) 


energy g of formation of chemical complexes values 
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which are within 20% of each other. Also, in each 
case, Wio=V>2- However, the values of W53 and 013 
differ in some cases substantially from the corres- 
ponding V53 and Vi3° This: isi probably rePated) to ‘the 
differences in the definitions of oe and Vi5 which 

are likeiy to be significant when the chemical complex 
Pe.anvoived, i.e. for i(i9)) = 2,3 vand 136) she - tracing 
of relations between w and v, however, must await their 
interpretation in terms of the basic interatomic 
interactions between the A and B atoms. We should 

also recall that the use of the expressions (5.7) and 
(5.10) for G' entails the assumption that the interac- 
tions oe and Vv; are small quantities - for binary 
regular solutions one usually considers the expression 
corresponding to (5.7) to be approximately valid for 
Jw, >| <2RT (see. $344)\2 ain Table 5.1, .some or | 53 | 
etc. are as large as 4RT. It is possible that these 
large values of Was (or vig needed to fit the data are 
partly reflecting the inadequacy of the approximations 
in which the combinatorial terms in (5.7) and (5.10) 
are treated. In this connection it is interesting to 
note that for the Mg-Bi system V53 and Vi3 are consi- 
derably smaller than Wo 3 and W137 SO that for this 
system the approximation (b) is to be preferred over 


the approximation (a). This is not surprising in view 
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of the discussion of the two approximations given 85.2 
and §5.4 - in Mg-Bi there is a strong tendency to form 
chemical complexes and utv (=5) is relatively large. 
In conclusion, it has been shown in this 
chapter, by four examples, that for some compound 
forming binary molten ayer ens the concentration depen- 
dence of the various thermodynamic quantities may be 
interpreted by assuming the formation of appropriate 
chemical complexes. The model then provides informa- 
tion on the free energy of formation of the chemical 
complexes Sale and on the interaction energies between 
the three species A, B and iy In wCNemSOLUCcITONe fic 
“will be realised that the model is a phenomenological 
one in which the concentration n3 of the chemical 
complexes occurs as an order parameter. It is not 
possible to explicitly demonstrate the existence of 
these chemical complexes from the thermodynamic data 
or from the zero wave number (g>70) “Limit. of the 
structure factors Sac (a) etc. However, if the chemical 
complexes are actually formed, then their equilibrium 
number and the geometric arrangements of the atoms in 
them should, in principle, be determinable from the 
knowledge of the structure factors at higher q values, 
just as the S(q) at higher q values reflects the mole- 
cular nature of a polyatomic fluid [Egelstaff et. al. 
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Finally, it is worth noting that the assump- 
tion of the formation of chemical complexes underlies, 
also, some of the recent discussions of the concentra- 
tion dependence of the electron transport properties 
of Mg-Bi and Tl-Te systems [Enderby and Simmons (1969), 
Enderby andCollings. (1970)), Cutler (28y2)yo stIn*such 
piacutaiene: of course, one has to make additional 
assumptions regarding the behaviour of the valence 
electrons on formation of the chemical complexes. In 
as much as the present work describes the concentration 
dependence of the thermodynamic properties quantitatively 
and, additionally, allows one to infer the number of 
chemical complexes at any given temperature and pressure 
from the thermodynamic data, it may be helpful towards 
a more quantitative interpretation of the transport 


properties of these systems. 
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CHAPTER 6 


THE STRUCTURE FACTORS FOR ARBITRARY q INTERACTING 


BINARY MIXTURES 


6.1 Introduction 


The work of the last three chapters shows that 
the concentration dependence of the long wavelength 
(wave number q>0) limit of the structure factors 
exhibit widely varying features depending on the type 
of molten system under consideration. In particular, 
we saw that the deviation of the concentration fluctua- 
tions Sac (9) from the ideal value (c)c,) reflects rather 
sensitively the nature of the interatomic interactions 
in the mixture. On the other hand, as discussed in §2.2, 
the sum rule obeyed by a structure factor and its short 
wavelength (q>«) limit are the same for every system - 
for example, aa td) and Sac (a) respectively obey the 
sum rules (2.22) and (2.21) and ei and Se 
C}Co- Thus we expect that the structure factors at 
intermediate gq also depend, in general, on the inter- 
atomic interactions in the mixture. We may recall, 
parenthetically, that the knowledge of the structure 
factors at non-zero q is crucial, not only for under- 
standing the usual X-ray and neutron scattering 
experiments in a mixture, but also calculating 


the electron transport properties of alloys of metals. 
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The calculation of the structure factors for 
arbitrary q is far from easy and, in fact, tne only 
calculations available to date seem to be for a me 
ture of hard spheres using Percus-Yevick approximation 
[Lebowitz (1964), Ashcroft and Langreth (1967), Enderby 
and North (1968)]. Here, as the name implies, one 


regards the system as a random mixture of hard spheres 


of appropriate atomic (molecular) diameters, so that the 


distribution of the two types of atoms (molecules) is 
entirely governed by their sizes. The model clearly 
corresponds to an athermal type mixture of Chapter 3 


ij-= 0 in Eqns. (3.25) “and (r4oyT It thus cannot be 


t It is to be noted that even within the framework of 


' the model, the derivation of the expression (3.42) 
for Gu and hence (3.44) for Sac (9), on the one hand, 


and of the expressions for the various structure 


factors S_-(q) etc. on the other, involve other mathem-— 


ce 
atical approximations. Hence, for example, the value 


of lim Sac (qd) from the latter calculations would, in 
q>0 . 
general, not agree with the expression (3.44) [with 


(oO) 


gees 0). *for'S..(0) a. TE, however, 0,/55= (Vy /V5 


(6 @: 
both give Sac (0) =e(1-e). For this case Sac (a) = Cii=G) 


forall gq and Suc (a) = 0, and we have a substitutional 


alloy at all q. 
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expected to give agreement with experiment at q+0O limit 
even for the simplest of real systems. For example, for 
3 (a) ence) 
Na- att Ze ~ 
Bohra ILOy 5 wis tla (5,/0y,) (Vi. /Vua yr 2 (Op, Ona! 
diameters of K and Na atoms respectively), the hard 


sphere model. gives, at c= 1/2, S_..(0) =024, whereas the 


cc 

experimental value, at 100°C, is Sac (0) = 0.55 (see Fig.3.4). 
The purpose of the present chapter is to present 

an approximate modification of the hard sphere treatment 


SsOnthat the q~ 0; limit of oo. “(G) eS omoOre aneaceor daiwa 


CC 
that for a weakly interacting system - rather than for 

an athermal or non-interacting mixture - and then to 
examine, by an example, the behaviour of the structure 
factors at non-zero gq on the basis of this approximation. 
The basic equations are collected together in 86.2 and 
the modification of the hard sphere model is described in 


§6.3. Section 6.4 gives the numerical example and is 


folgowed bysa brief descussion of the results andveners 


significance (86.5). 


642  Basic*Equations 


622.4 “Direct Correlation Functions sancrotruccurcsnactons 


There are essentially two approaches to the 
calculation of the pair correlation function g(r) of 
Aebiguid sor sot Gig it) for a Mixture) ihe vous, 
developed variously by Bogoliubov, Born and Green, Yvon, 
and Kirkwood [for a review and references see, for exam- 


ple, Cole (1967)] results in*a heirarchy of differentio- 
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integral equations which express the pair correlation 
function in terms of the interatomic potentials and 


(3) 


the peated correlation functions g , tiercriplet 

in terms of quadruplet and so on. In actual calcula- 
tions one truncates these BBGYK equations at the first 
by making some convenient approximation for the triplet 


(3) 


correlation functions g such as the super-position 
or convolution approximations. The second approach 
involves instead the substitution of an approximate 
expression for the direct correlation function, first 
introduced by Ornstein and Zernike. Of several appro- 
ximations available, the Percus-Yevick (PY) theory 
{[Percus and Yevick (1958)] has proved most convenient 
and we base our calculations Seweten 

For a binary system, the direct correlation 
functions Sia? = Saas: i,j = 1,2, are defined by 


the relations 


2 
, Sera 
his (r)=C,.(r) +n wey cy [Bag (2 Cos ([x ec Rote (onth) 


where n is the number density, Cy is the concentration 


of atoms of type & (2 = 1,2) and 
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Denoting the Fourier transform of hy, tr) by his 


his (q) = ine (r)etd-t a> ; (6.3) 


and similarly the Fourier transform of C,4,) by C5 (ale 


one obtains from (6.1) 
2 
hj 5 (q) = Ci; (q) + Be coh; (q)C,, (a) (6.4) 
Remembering the definition of the Faber-Ziman 


structure factors a... (q) ;4Eqn..(2.1), wonehas 


1)j 
iq-r 3 
an (a) = 1+ nf tg,5 (2) - lje™4 Geer 
=] + nh, 4{q) A (6.5) 


etl from (6.4) and (675), 
a1, (a) = Do *(q) [14ne, (C,, (a) -C5, (a)) 
Fen cS (GC) Squceeuay ei 
seta) — Di (a sila (Ce layla nia) 
(6.6) 


2 
SntceKcy, 1g)¢,, (Gg) C7, (a) 1 


a. = pi {qi} [1 = n(e, Cy, (ah t'e5C5, (a) — C5 (a)) 


ae 


+ nics GCppahess ia) eal} qua) 


where 


2 2 
D(q) = l-n (CC, 4 (a) tenCa9 (q)}tn™c) Cp (C1 (a) C55 (q)-Cj 5 ep) yee 


(6.7) 
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For a single component system the Eqns. (6.1) (6.4) and 


(6.6), of course, take the well known simple forms 


h(r) = ctr) + afn(etye( [rer [yae? 
h(q) = C(q) + nh(q)C(q) 
and 
a (q) = S (q) = TREY ° (6.8) 


The expressions for the number-concentration 


structure factors Sy (2) etc. for the binary mixture in 


terms of Ae eeu are readily written by substituting (6.6) 


into (2.17). In the following it will be convenient to 
work with the concentration fluctuation structure factor 


Sac (a) and the combinations 


2 
6 (q) = Suny (a) - [Sue (4) /Sa¢ (a) J Sac (a) (6.9) 
and 


A(q) = -Syic (4) /Sa¢ (4) : (6.10) 


We observe from (2.32) and (2.33) that 6(q) and A(q) 


are defined such that 


640) = 8 = AY) kT kp 


and A(0) = 6, the dilatation factor. Substituting (6.6) 


into (2.17), one obtains 


iL 
Bio) OS ere an es a a ee C6228) 
1=n (62, , (a) +5099 (a) $20,650 9 (a) ) 
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peal ge me ta) Lo (Cy aye) ta) ) 65 (C.) (aye tao (6.12) 
and 

2 
€1¢y[1-n(c} Cy, (q) +e5C5, (q)+2c,¢,C, 5 (a) 


aefes eUlas 2 2 
I=n(cjC, 4 (q)+e5C,5 (q)) +n C1 CSo (C14 (4) C59 (Gq) -Cjy 5 (a) ) 


(6.13) 


On using (6.11) and (6.12), Sac (a) can also be written in 
the form 


ae 


Sac (a) = q 
pai cs Meer ate Whe Fe pik iets we CEL aeateo T(oy 


which will be found useful later. 


6.2.2 P-Y Approximation and Hard Sphere Model 


In the P-Y approximation for a mixture, the 
interatomic potentials $5, (2) are related to Sea by 
the relations [Percus (1962), Lebowitz (1964) ] 

wl Sool Gy -B6o.. (r) 
i) = = 1) 
93; (x) le 1) eae eels (62.05) 
where £ = 1/k,T. Equations (6.15) together with (6.1) 


and (6.2) are sufficient to determine eee or 954 (x) 


a2 De ks are known. 


For a mixture of hard spheres of diameters Oy 


and G5, one has 
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which imply, from (6.15), that 


o2...\(r):1=..0 Lousera< ake) 
and 


C..(r) = 0 LOr, = 2 > aioe.) P (6.17) 


With (6.17), the solution of the Eqn. (6.1) for Shee 
is [Lebowitz (1964)] (x =r-)) 


3 : 
a,+ b.r+dr Ge APO sl we le 2 


Rey) 


-C,,(r) =a 


12 


a,+ [bate Waa as 7 ae ete (0, +05) 


(6% 18) 
where we have assumed, for definiteness, that 0, 2 9) 
and where ass bi, a,b, d are constants which are func- 
tions of O11 Ips the number density n of the atoms and 
the concentration of two types of atoms. For brevity 
they will not be written here explicitly. When 0, = 55 
the various constants are independent of concentration 


; t 
and a, = apr by = bo = b, so that 


Cy!) = Co (r) = Ci2(*), ror 0, =% - 
2 ee es Sen ere ranean he i Chak eB eels A os 
t As to be expected, each Coy) is, ‘then, juste Gof a 


one component hard sphere fluid, the expression for 


which was first obtained by Wertheim (1963) and Thiele 
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Remembering the relations between the Fourier transforms 
e532? of sie and the various structure factors, 


one sees that for hard sphere mixture of equal diameters 
Sac (a) = C]Co , A(q) = 0 


and aj3 (4) = @5 5 (q) = Ay (a) 


so that one has a simple substitutional mixture. 

The calculation cy ta) from (6.18) for any OF 
and Sor involves integrals of the type [sin qr rvdr 
with n = 0,2,...,4, which can be evaluated in closed 
form. The resulting expressions for C45 a) and hence 
the structure factors have been obtained by Ashcroft 
and Langreth (1967) and Enderby and North (1968). ee 


~ they are lengthy we shall not reproduce them here. 


6.3 A Modification of the Hard Sphere Model 


The interatomic potential consists of a strong 
short range repulsive interaction and a relatively long 
range weak attractive interaction (see Fig. 6.1). The 
hard sphere model takes into account the former in 
the form of an infinitely strong repulsive core but 


neglects the attractive interaction. An approximate 
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method" 


of taking the latter into account is suggested 

by the following argument. | 
We notice from (6.15) that if $54 (°) = 0, at 

a given r, then provided G54 (2) is non singular, 

C,4 (0) =O... Further if C,5'°) are zero for all r, then 

eqn. (6.1) gives cha = 1 for all r [ideal gas or 

hard sphere model of zero diameter!]. This suggests 

chats is [8o; (x) | << 1, we can set, approximately, in 

(6225) 9447) ~1l and thus obtain 


fate | . 
Cy5(r) = 8655 (r), lds.) | << 3g . (6.19) 


In analogy with the case of one component system, this 


approximation may be referred to as the random phase 


oy 


The method is an extension of that proposed by 
Woodhead-Galloway et.al. (1968) for a pure fluid. 
While this work was in progress it was pointed out 

to us by Professor March that a similar extension has 
been suggested by Woodhead-Galloway and Gaskell (1968). 
Our procedure of specifying the repulsive core and the 
attractive potential (see below) differs however from 
that used in these works and the discussion of Sac (a) 


etc. presented here has not been given previously. 
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approximation (RPA) for the mixture. ! Because 9; (x) 

is not small for small interatomic separation, the 

RPA is actually quite a poor approximation for liquid 
media. To take into account both the strong short range 
repulsive interaction and the weak long range attractive 


interaction, we assume that C;4 2) has the form 


~ CHS us 


where Cyy te) are the direct correlation functions for 
the mixture of hard spheres of appropriate (see below) 


diameters oF and 5 and 


t For a pure fluid, the Fourier transform of the exact 


Born-Green-Yvon equation connecting g(r) to $(r) and 
the triplet .correlation function is) (q (750) a isee eng. 


Brout (1965) ] 


B k.g P . 

< eS GN ie) Oe iO eae (6.19) 

PgP-g N D k kPg?-K-q 

N  ig-R 
where acie=ait), ble a Ry being the position of the ee 
gq tel 

atom. The RPA corresponds to neglecting all terms in 
the sum in (6.19)' except k = -q, whence 
S(q) = i. Oe nt hea ene ee C69) 

NG gd q 


which implies, comparing (6. 8)2,.(C(Qyt = -8n$,- 
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where ri is the position at which acs) is minimum. 


The approximation (6.21) implies that beyond re 


wy the 
13 


potential can be considered to be weak, since use of 


Ges 2 1) pedslivensOne-s Ounces 


das 28) 


as in RPA. 

The choice of a non-zero Xqg (tl, LOU aa ri 
needs some comment. First, if eee) is set equal to 
Becton ror 7 <riye then this abrupt cut off gives rise to 
spurious harmonics in C54‘) and hence in various struc- 
ture factors. Secondly, as discussed by Gaskell (1970) 
FOGeagpure £luid,, a nNon=Zexro Kqa'®) Ln thesecore reqvon 
is necessary to satisfy the physical requirements that 
thespair .cistrisueion Euneciens aie) (calculated from 
6.20-6.22) be zero for pees ANG Boa! To fend esuch sa 
aye a difficult mathematical problem and Gaskell 
in his work on argon also took y(r) to be a constant, 
Say Xo» inside the core region. He adjusted the value 
of Xo $0 that the average of g(t) “over the cores =zera, 
ret oe i g(r)dr = 0. This gives for argon at 86.5°K, 


— 
BX, Ei =1.67 compared to the value ot Ses = -1.39 from 


(6.22). Our suggestion for y¥(r) has, perhaps, the 
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advantage that it requires no additional work to deter- 
mine Xo and that it four y(r)] is smooth over the whole 
range of r. We note parenthetically that in Gaskell's 


Work, x(t) =?¢(r) for’ rr > ry and x(r) = Xp for rs ro 


where ro (= o) is the value of r at which g(r) = 0. 


6.4 Numerical Calculations 


6.4.1 Calculations for Liquid Argon 


Before presenting calculations for a mixture, 
it is useful to consider the application of the above 
approximation to a pure liquid. We choose liquid argon, 
for which the interatomic potential is usually given in 
the Lennard-Jones form: 

5 eerene OF Tae 
g(x) = 4e/(2)  - (>) d (6.24) 


For argon, 


Soon ee ee ee (6225) 
k oO 

B 
We make the calculation at 85°K at which molar volume 
is 28.34 cm and n = Deda 0" + fem. The functions 

: : : 1/6 

@(r) and x(r) are depicted in Fig. 6.1 (x =2 ‘ re 
One may physically expect the appropriate hard 
core diameter o to be somewhat less than the value 


(r_) at which the potential is zero. Although 
fe) 


the fit for the structure factor, S(q) (=6 (q)) 
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Pig.'6.1.. Sketch of the model potential showing 
¢(r) of Lennard-Jones form and x(r) 
as well as the hard core. 
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with experiment was not unreasonable for the choice 
o= 2X5 a somewhat better agreement with experiment is 
obtained for the choice o = 0.98 lo: The theoretical 
and experimental results are shown in Fig. 6.2. We 


observe that the agreement with experiment is not 


unsatisfactory. 


6.4.2 Mixture 


To examine the effects of the long range attrac- 


tive interactions on the various structure factors of a 
binary mixture, we shall here consider, for simplicity, 
a mixture in which $4 (Fe $55 (r) and 42 (¥) have the 
form (6.24) with the same ror so that one may take, to 
a first approximation’, the hard core’ diameteusr for) the 
two species to be equal (o, = NGO say). Further, 
the three Dee hee are then proportional to one another, 
and similarly the three eee and hence their Fourier 
transforms Ges Denoting X41 '@) = x(q) and intro- 


ducing two dimensionless parameters y and: 


| E12 E11* €o27 2£12 
a : 
mate i 


’ (6.26) 


to measure relative differences in the (minimum) depths 


of o55(r), one has 
a 


: This may be expected to be reasonable if the diffe- 


i i see below). 
rences in various €i5 are small ( 
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ql 


Figs (6:2iqModel~$iq),. (solid line) for </kg= 1230Ky 
Yo= 3.4x10°° cm and o = 0.98 ro compared 
with experimental S(q) for Ar at 85°K. 
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X71 (G) = x(q) , Xy2 (a) = (1- y) x(q) 


and Xoo 'G) ots (22 y tly tq). (6.27) 
Finally, recalling (6.20), one has 


Ci4{a) = Suga) - 8x, , (a) (6.28) 


where Cig (a) is the Fourier transform of the direct 
correlation function Cig §¥) for a fluid of hard spheres 
of diameter o. 

Substituting (6.27) and (6.28) into the 
expressions (6.11), (6.12) and (6.14) for 8(q), A(q) 


and S..(q), one obtains 


ue 


@(q) = —— | (6.29) 

1-nC,,, (q) + nBx(q) [1 - 2ycg + ARco] 
A(q) = -n@(q)xlq) Bly - 2Rc 4) (6.30) 
and 
i c,c 
S.,(q) = x2 4 (6.31) 
oo A“ (q) 

1+ 6(2nx(q)& )cjc, - CC costa 

i ia 2 Enon then, from (6.26), y = and the expressions 
for 0(q) and A(q) take the form (y = ®) 
8(q) = Sig ese ee ie ee {6.32) 


1-nC,,. (q) 4 ney iq) (l= 2Rc, Cc.) 


A(q) = -n@(q)Bx(q)A(1- 2c5) .- (6.33) 
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We note that this case corresponds to systems where 

the volumes Cae and ey) of the two pure species 

are equal. Indeed, the expression (6.33) for A(q) 

has the same concentration dependence as the expression 

(3.31) for 6 © A(0)) for a conformal SOLUELON ,~ bream 

the latter one sets eae oa: Finally, for a mixture 

which satisfies the substitutional bypothesis, £11 = £92 

E,2 and hence y=R=0, so that 6(q) is independent of 

concentration, A(q) =0 and Sac (a) =] eo, as expected. 
We observe from (6.33) that at c = %, A(q)= 0, 


so that, for this case (Q= y, Co = 3), the value of 


Soc (a) depends on the single parameter 


nx(q) >A = ot (6.34) 


which defines w(q). If we use the interatomic potential 


fori liquid. argon to-calculates vq) ,we find etome, (0) = 


x (q=0) 
co 

n n 2 
(ap) x (0) = 4m (ea) | reGe) ie isla 

B B 0 

32 30P on a ero ee 
ae tee ore. a = T v (6.35) 

and hence w(0)/k,T = -1624 @/T. For mixture having this 


value of w(0), Sac (0) would become infinitely large, 
T= 85°K if @ = 0.105. For &> 0.105, S,,(0) would be 


negative which is not allowed since Sac is non-negative 
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by definition. In other words, for @ = 0.105, the 
Gritical temperature of mixing is 85°K. Now it may be 
seen from eqns. (6.30) and (6.31) that if |@| and 

ly| << 1 (or, of course, if Percerediya, aces aU) shee ie 
the differences in the pairwise interaction energies 
are small, then since 6(0) © 0.06,-the term involving 
iin Che denominator of (6.31) can be ignored at all 


Concentrations for calculating S..(0) wand hence: (o-o1) 


CE 
gives 
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kp deer 


which is just the expression for Sac (9) derived in 
Chapter 3 from the expression for the free energy of 

a regular solution (in the zeroth approximation) or a 
conformal solution. As remarked earlier (see p. 127) 
we cannot in general expect that the approximate 
methods: for calculating the Sac qd) etc, £oreanrbitrary 
g, would, give for Sac (9) etc., the same expressions as 
derived from the appropriate statistical mechanical 
theories of mixtures described in Chapter 3. It is, 
therefore, not unencouraging that the method suggested 
here gives, at least for the case vacua! el ewrang 
iGimand ind <<ilroan expression for Sac (a) whachein the 


limit q+0 agrees with that derived from statistical 


mechanics for this case. 
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In Figs. 6.3-6.6, we present some results of 
numerical calculations to illustrate the variation of 
the different structure factors with wave number q. 
For definiteness, we took throughout for y(r) and o 
the values appropriate to liquid argon as given in 
§6.4.1 and assumed for the mixture that @= y. 

The results in Figs. 6.3 (a and b) are for 
(UY = 0.10 and at T = 85°K which (for this value of f%) 
as discussed above is close to critical solution tem- 
perature of the mixture, so that Sac (0) is very large 


as q>0 for c)=c,=%. Fig. 6.3a shows 6(q), Sac (a) 


and A(q) for Coe and Co= 0. 2 [A(q), in our examples, 


is zero at oie ee and is in general quite small at other 


concentrations since we are considering here only the 
mixtures with CRU iy Sera Fig..6.3b shows Sg for 


the same two concentrations. 


Fig. 6.4 gives similarly the results at an 
elevated temperature T= 5 ae using the same & namely 
Gearon. 

As will be clear from the definition of &, 

a positive value of @ implies that like atom pairs 
have lower energy than the unlike atom pairs. To 
illustrate the behaviour of the structure factors 
when unlike atom pairs are energetically preferred, 
Fig. 6.5 gives the model structure factors. for 


@ = -0.2 at T= 85°K. As expected, S,.(q) for q>0, 
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Fig. 6.3a. Model 6(q), Scc(q) and A(q) (inset) for two 
concentrations (c = 1/2, c=1/5) taking 
R= 0.10, T = 85° (i.e. critical parameters). 
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Fig. 6.4a. Model 6(q), Scc(q) and A(q) (inset) 
at an elevated temperature: T= 3/2 T 
with R= 0.1 (0 = 85"). = 
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Fig. 6.4b. Model aj3 (q) for the same parameters as 
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Fig. 6.5a. Model @(q), Scc(q) and A(q) (inset) for 
| R = -0.2, T = 85°K for two concentrations. 
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Big...6.5b. Model va: 
Fig. 6 


for the same parameters as 
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is now less than the ideal value (c) C5) in contrast to 

the case fR>0. Also, it is found that the oscillations 
in Sac (4) damp less rapidly when &< 0 than when A> 0 

(Fig. 6.6). 


6.5 Discussion and Concluding Remarks 


The major feature of the results depicted in 
Figs. 6.3-6.5 is that the effect of the difference in 
the interatomic potentials on the structure factors is 
primarily confined to the values of q in the range 
Oo <g AS z dor where qo is the position of the main 
peak in 0(q) or a; ; (a). For example, for q less than 
about 5 doe Sac (a) is quite different when #%< 0 from 
that when &{> 0, while beyond this q, Sac (a) =C1Cor 


irrespective of the value R. The partial structure 


3 


another (except that A471 (4) = ago (a) for Cy =%, since 


factors a4) for tee Gor are different from one 


in our example A(q) =0 at this concentration) and depend 
on the concentration, the dependence of ais on both q 
and concentration being quite different when @< 0 and 


2 
when @&> 0. In contrast, for q 2 5 Gyr one has appro- 


ximately 11 (4) = ag2 (a) = ay (a) and the dependence on 


concentration is quite small. In other words, the a;5 (a) 


here approximately obey the substitutional hypothesis. 
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Fags! 6.6. Sac (4) for’ two values of Rat T= 85°K. 
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In the calculations presented in Figs. 6.3- 
6.5, we have concer d the mixtures with Wee Macd 
so that. the hard core diameters Oo; for the two types 
of atoms are equal. When the two diameters are not 
equal, one would still expect that the effects of the 
(weak) long range interactions are confined to q-values 
less than a value qo where Sate ae For de ae the 
structure factors may be expected to be approximately 
those given by the hard sphere approximation with 
different Oj 

Our work thus shows that for weakly interacting 
mixtures - that is, where the differences in the inter- 
‘atomic potentials are small - the structure factors in 
Ghatereanges0 < qi< q, can be significantly different 
from those calculated from the hard sphere model, the 
difference being dependent on the relative strengths 
(magnitude and sign of #¥ and y) of the potentials. 
Paehoagh a comparison with experiment of the structure 
factors for arbitrary q is at present not possible 
because of lack of experimental scattering data, the 
conclusion arrived at here are nonetheless of consi- 
derable interest. 

For example, in the theory of Faber and Ziman 
(1965) of the electrical resistivity of mixtures of 
simple metals, one has essentially to integrate the 


scattering function I(q), Eqn. (2.27), with a weighting 
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factor aa) from q = 0 to 2dpr where dp is the Fermi 

wave number of the electrons. The calculations have 
hitherto been performed either using the substitutional 
hypothesis or the hard sphere structure factors with 
different o. For multi-valent metals, 2p >q, so 

that the integration extends beyond the first peak 

and hence the use of hard sphere structure factors would 
form a good approximation. However, for monovalent 
metals 24p< dos so that the modifications of the hard 
sphere structure factors at low q would, in general, 
have to be taken into account to obtain fully quantita- 
tive results. Since the interatomic potentials in metal 
alloys are not so well known, a useful procedure might 

be to choose a simple analytic form for X44 (F) and fix 
the parameters such that the long wavelength limit of 

all three structure factors agrees with those derived 
from thermodynamic data. The band structure and Madelung 
energy Ae depend on the structure factors and a similar 
procedure might be advantageous there too. 

Finally we recall from Chapters 4 and 5 that 
the long wavelength limit of Sac and other structure 
factors for compound-forming molten systems vary with 
concentration in a spectacularly different manner than 
for the conformal or weakly interacting mixtures. In 


these (compound-forming) systems the associative tendency 


y an 
Ae — 
ae ' , 
ue , yk 
can 4 . ae } rial Y : nm 
\ n oa ie i ne as 
y 0 as a 
" re ) (0) a ee ae eye ; 4 
1 hy | Sieh ian hy Tw ha My » re : ow J £ Nowe a 
uf TA he af Le ee ms Fi 
— an i ta y i iy a! } = rit mY ay) 4 
' Pie eta ' ‘ vi Ae vi 
shi ; 7 Di he eae lll “ 4 a on Ps ( + i a ‘ 
en PR aoc a4 Lae: “ mio! 0 
P Ue my ie 
..) 
bed hohe. en tyre ot ae ad aes 
». ‘bale webheud oi iw aa a we 0 a ee ae” a -" ¥ 
wea ce oe hea scent a Oe De iy ram et 4 i! iy DHA be 8 Ra he <a es ? he 
i ] fs Pa | ; 7 | j om \ | oe an | ie 
val re ee cn Cx — dire 7 we oh f er ek aati 
ut ' | eh a re Hei: * ota peng HLT te veep a oe y yet et Le 4 ret me hg siti, a 
gy eee ae a : 
f nt Lo". iad Se ek eb ey O14 te ae “i ana 
mete mag 
he, 
i r ' 
o% » & ( * wie e ” , 1 
time eo LASS, B War Ste 2 Oly 
i 7 y aL pre) ; 1") ay dt ra 
ayo, “OR . ew » ert oad 293g qe 
j wae i f 
i ie f & aie | a wy ; 
‘ \ eer eee a: ee ee &"wW CC. : phys ee 
tf Bo 9O BHohieo tebe oh? cade Oy tan ° ai 
[ } i 
i i i‘ ard ; x Cm) ar 
ry i ata im 
hy 
ee TT a Low oe Wee ee 5 i ae pend ay 
mM PA 5° Pe an ’ x 
£ nN i i ~ 4 : * : + é neh 
ee Roa iby d i ie , va. he a wns AN. o- a ae 
J Ai! Gad tne: ee hie NS tee MT, ? yh! Se a re Agee © nite or, ~ hal ID 
al i Mi LNA \ ae i 
t nt iad % en in pe - ba 7 aan Can i Fr Wes 
be weg DE bah rie A Sy ead ia ab te tee) eee © 
, i 1 / i Me kaa fai ie mh 
, ee ee eee bs 
nh | Ree SPs: a sh Ta rent tes ay, 
f %, ~~ ri 7¢ i. t 
H i ; i 
r a re J 
y é Ce fet ‘ hy 1k orcs iF twos ‘os! oa ts de 
i. n } ay i) io ay 
i i v F ai tS 
hy "| P| ns h 
Sees iPr e. ye viteh, tt ine i , Ho ries ie tes ; es yah @ 
i i iy 
aa a 
; Te) Pe; Wateot Bi , 
tne 
i rs e Wea " 
“ hoi id hhh ey i 
Pn 
iF fa 
ae 
eee a 
eee OTR 
iy 
ape | Oh 
Thiel beat le 
hiv fare eae 
a ol ls WY ol ; “ er 
és ie Pu ¥ 
: os Wy oF ‘ ; 
x ae ob aay LA ed dtAt 
TS aa . 
a ' vy ' 7 rar a 
' a i c ahi ai 
i) me: i’ : siete ge ct nd re 
ft ' ~ 
De 
, i 


jdt Sh Come any 
i oie a Yi 


ESG 


(characterised by stoichiometry) between the atoms is 
very strong and the metiod of calculating Sac (a). etc. 
presented here is not sufficient. Recently Ruppersberg 
and Eggar (1975) have measured from neutron scattering 
Sac (@) for LiPb mixture at the compound forming concen- 
tration Li,Pb, and Sac (a) shows oscillations which are 
of much larger amplitude and which damp less rapidly 
than Sac (a) of Fig. 6.6. It was shown in Chapters 4 and 
5, that the concentration fluctuations Sac (9) and other 
thermodynamic properties of the compound forming molten 
systems may be quantitatively explained on the basis of 
POCHAt ion of appropriate chemical complexes. Recent 
work [Bhatia and Ratti, to be published] indicates that 
the idea of complex formation can also satisfactorily 


explain the neutron scattering data. 
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APPENDIX 1 


The formulation of Chapter 4 does not work 
with G', introduced in Chapter 5, but rather with the 


derived quantities Ne which are related to G! by 


AO) ae 
RT 2n Vs SoG ce auar Be G; RT gn (n, /n) 


(aG'/on;) - RT gn (n, /n) : (Az D) 


prop iie 


Qn Y,; are zero if the ternary solution of A, B and Baty 
is ideal. In (A.1) the suffix n' means that in differ- 
entiating with respect to nis all other n, are kept 
Constant. Because of Che extensive property of Gy 


(aG'/on, ) and hence {£n y,; satisfy Gibbs-Duhem type of 


relations 

. n, (£n REELS yee ae = O05. (AR) 
Since G’ = ) re OG EE , one thas’ from °(572) 
G, = —nag t+ RE ) n, gn (n,/n) + RT ) n, &ny; (R38) 
and the equilibrium condition (5.5) gives 
(ny y,/n) # (nyy/n) © = (n,y3/n)expl-g/RT] , (A.4) 


(A.3) and (A.4) are the same as given in Chapter 4. 


The advantage of using uA is that like (A.4), some 
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of the formulas can be derived with a minimum of 
algebra and expressed compactly. For example, the 
activity an is simply (y)n,/n). However, since in 
theoretical work it is the G' which is first calculated 
and Y, are the derived quantities, we have not used 
here the formulation of Chapter 4. The two treatments 
are, of course, equivalent. 

For the approximation (a) for G', expression 
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